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Abstract 



The flavor singlet axial charge has been a source of study in the last years due to 
its relation to the so called Proton Spin Problem. The relevant flavor singlet axial 
current is anomalous, i.e., its divergence contains a piece which is the celebrated 
6^(1) anomaly. This anomaly is intimately associated with the rj' meson, which 
gets its mass from it. When the gauge degrees of freedom of QCD are confined 
within a volume as is presently understood, the Ua(X) anomaly is known to in- 
duce color anomaly leading to "leakage" of the color out of the confined volume 
(or bag). For consistency of the theory, this anomaly should be cancelled by a 
boundary term. This "color boundary term" inherits part or most of the dynam- 
ics of the volume (i.e., QCD). In this thesis, we exploit this mapping of the volume 
to the surface via the color boundary condition to perform a complete analysis 
of the flavor singlet axial charge in the chiral bag model using the Cheshire Cat 
Principle. This enables us to obtain the hitherto missing piece in the axial charge 
associated with the gluon Casimir effect. The result is that the flavor singlet 
axial charge is small independent of the confinement (bag) size ranging from the 
skyrmion picture to the MIT bag picture, thereby confirming the (albeit approx- 
imate) Cheshire Cat phenomenon. 
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Chapter 1 
Introduction 



The constituent quark was proposed to explain the structure of the large num- 
ber of hadrons being discovered in the sixties [|TJ . Soon thereafter deep inelastic 
scattering of leptons off protons was explained in terms of point-like constituents 
named partons Q. The analysis of the data by means of sum rules led to the 
conclusion that there was an intimate relation between the partons and the ele- 
mentary quarks. Various models have been developed to understand the struc- 
tures of light hadrons and their interactions in terms of quarks ||. They were 
built on the basis of low-energy hadron phenomenologies, particularly, (i) approx- 
imate SU(3) flavor symmetry and its explicit breaking, (ii) Okubo-Zweig-Izuka 
(OZI) suppression rule for flavor changing processes JiJ, and (iii) chiral symmetry 
realization with its spontaneous breaking pattern. 

The birth of Quantum Chromodynamics (QCD) and the proof that it is 
asymptotically free set the framework for an understanding of deep inelastic phe- 
nomena beyond the parton model [[J. However the fact that QCD confines does 
not allow a solution of the theory in the strong coupling regime and therefore 
new models had to be developed to describe hadron structure which realized the 
phenomenological principles mentioned before but in a manner compatible with 
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the dynamical principles of the theory. This scheme of confined isolated quarks 
and gluons, has a strong relation to the valence quark model for hadrons || . The 
valence quark model was developed further to a non-relativistic quark model by 
De Rujula, Georgi and Glashow |?J initially and exploited phenomenologically by 
Isgur and Karl ||, and into a relativistic quark model framework by Chodos et 
al. and De Grand et al. known under the name of MIT bag model ||. 

Although the MIT bag model was successful in describing the properties of 
the nucleons, it was not pertinent due to lack of the chiral symmetry. In order to 
incorporate chiral symmetry into the model, the pseudoscalar mesons had to be 
introduced in this framework [[K|]. The resulting scheme, the so called chiral bag 
model, was constructed with meson fields which are restricted to be outside the 



bag 11 



Chiral symmetry, a property of QCD with massless quarks, has been instru- 
mental in the description of hadron phenomenology. So much so, that Skyrme 
realized its importance and wrote down, much before QCD, an effective theory 
for the strong interactions, in terms of pion fields only, describing a unified theory 
for baryons and mesons ||12|| . Only many years later his tremendous intuition was 
appreciated and his ideas justified from the point of view of QCD |13 . 



The chiral bag model incorporates in a unified description the statements 
above. It is defined by means of a QCD lagrangian inside the bag and by a Skyrme 
type theory outside, properly matched at the surface to preserve the classical and 
quantum symmetries. These formulation leads to a intriguing principle referred 
to as the Chesire Cat Principle (CCP) 0, The possibility of formulating a 
physical theory by means of equivalent field theories defined in terms of different 
field variables, leads to this construction principle for phenomenologically sensible 
and conceptually powerful models. This principle states, that physical observables 
obtained by means of equivalent theories defined in a certain space-time geometry, 
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adequately matched at the boundaries, are independent of the geometry. In 1+1 
dimensions fermionic theories are bosonizable jnj and the CCP can be made 
exact and transparent. In the real four- dimensional world, bosonization with a 
finite number of degrees of freedom is not exact. However based on the unproven 



"theorem" of Weinberg |17j , it seems possible to argue that the CCP should hold 
also in four dimensions, albeit approximately. 

Quantum Chromodynamics (QCD) is the theory of the hadronic phenomena 
[||]. At sufficiently low energies or long distances and for a large number of 
colors Nc, it can be described accurately by an effective field theory in terms of 
meson fields [jT3| , |18| . In this regime, the color fermionic description of the theory 
is extremely complex due to confinement. However the implementation of the 
CCP in a two phase scenario called the Chiral Bag Model (CBM) has proven 
surprisingly powerful |19| |. The CBM is defined by dividing space-time in two 
regions by a hypertube, that is, the evolving bag. In the interior of the tube, the 
dynamics is defined in terms of the microscopic QCD degrees of freedom, quarks 
and gluons. In the exterior, one assumes an equivalent dynamics in terms of 
meson fields, i.e., one that respects the symmetries of the original theory and the 



basic postulates of quantum field theory p7| . The two descriptions are matched 
by defining the appropriate boundary conditions which implement the symmetries 



and confinement |M|, [L9| . What this does effectively is to delegate all or part of 
the principal elements of the dynamics taking place inside (QCD) the bag to the 
boundary. We will see that this strategy works quite efficiently in the problem 
at hand. In this scenario the CCP states that the hadron physics should be 
approximately independent of the spatial size of the confinement region or the 
bag fl4| . This realization of the principle has been tested in many instances in 
hadronic physics with fair success |L5) . 

There is one case, however, where the realization of the CCP has not been as 
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successful as in the other cases, namely, the calculation of the flavor singlet axial 
charge (FSAC) of the nucleon. Indeed in the previous results |2IJ |2T], the 
CCP was realized only partially as it seemed to fail at certain points such as for 
zero bag radius. It is the leitmotiv of this work to remove this apparent failure. 

Experiments using polarized electrons on polarized targets were carried out at 
SLAC f23|. Further information came from the SLAC-Yale group with fascinat- 
ing implications about the internal structure of the proton p4| . More recently, 
the European Muon Collaboration (EMC) obtained very extraordinary results by 
the scattering of a polarized muon beam with energy 100-200 GeV on a longitu- 
dinally polarized hydrogen target at CERN All these results point towards 
a new scenario in hadronic structure dominated by a quantum anomaly. To be 
more specific, the unexpectedly small asymmetry found by EMC implies a strong 
violation of the so-called Ellis- Jaffe sum rule and therefore implies that the 
polarization of the proton is not carried exclusively by the valence quarks. This 
problem is called the Proton Spin Problem PTfl . 

The EMC result and this problem are now believed to be resolved through the 
beautiful relation between the flavor singlet axial charge and the axial anomaly 
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a°(Q 2 )=^-N F ^ 1 ^g(Q 2 ), (1-1) 

where a°(Q 2 ) is the flavor singlet axial charge measured by EMC, AE the quark 
polarization, and Ag(Q 2 ) the gluon polarization. 

The aim of this thesis is to show the full consistency of the CCP in the hadronic 

world for the case of the Proton Spin, which was not satisfactorily established in 

the previous results in this direction [^, [H], |22]jp|. 

1 Note that in these papers, they have shown that the CCP holds for non-zero bag radii but 
it failed when the bag radius shrank to a point, implying that in the model studied, the pure 
skyrmion and the MIT bag did not have the equivalent structure required by the CCP. 
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In the CBM, the scenario of how the CCP is realized - which is the central 
issue of this thesis - is very intricate. As stated, the flavor singlet axial current is 
associated with the anomaly and effectively with the rf meson. Thus, besides the 
pion field of the conventional effective theories which accounts for spontaneously 
broken chiral symmetry, the correct treatment of the flavor singlet axial charge 
requires minimally the inclusion of a field describing the rf meson. 

The intricacies of the hedgehog configuration and its relevance to the fraction- 
ation of baryon charge and other observables have been extensively discussed |3(| 



and fairly well understood |3l], [32|]. They will be implemented in the present 
calculation without much details. Moreover the inclusion of the rf meson car- 
ries subtleties of its own. The vacuum fluctuations inside the bag, that induce 
the baryon number leakage into the skyrmion [BO], also induce a color leakage if 



a coupling to a pseudoscalar isoscalar field is allowed f33j. This leakage would 
break color gauge invariance and confinement in the model unless it is cancelled. 



As suggested in [33], this color leakage can be prevented by introducing into the 



CBM Lagrangian a counter term of the form 
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C CT = j> d(3 i^n^Trlnf/t - TrlnC/) (1.2) 

where Np is the number of flavors (here taken to be =3), (3 is a point on a surface 
E, n M is the outward normal to the bag surface, U is the U(Np) matrix- valued field 
written as U = e ln ^e lv '^° and the properly regularized Chern-Simons current 
K" = & vaf3 {G a v G a afi - lg s f abc G a u G b a G c p ) given in terms of the color gauge field GJ. 
Note that the counter term ( |1 . 2| ) manifestly breaks color gauge invariance (both 
large and small, the latter due to the bag), so the action of the chiral bag model 
with this term is not gauge invariant at the classical level but as shown in |33| . 
when quantum fluctuations are calculated, there appears an induced anomaly 
term on the surface which exactly cancels this term. Thus gauge invariance is 



restored at the quantum level. 

The equations of motion for the gluon and quark fields inside and the r( field 



outside are the same as in [EH, |2l||. However the boundary conditions on the 



surface with the inclusion of eq. (1^2) read 

n-E a = -^n-BV (1.3) 

87T7 

n x B a = ^i^h x EV (1.4) 

and 

in • (^775^) = fn ■ d V ' + ^fh ■ K (1.5) 

where E a and B a are, respectively, the color electric and color magnetic fields. 
Here ip is the QCD quark field. 

The full Casimir calculation of the gluon modes, which is highly subtle due to 
the p-wave structure of the rj -field, has to be performed to get the CCP for the 
flavor singlet axial charge. Here we would like to side-step this technically difficult 
procedure by first assuming the CCP in evaluating the Casimir contribution with 
the color boundary conditions (|1.3j), flTjj ) and fliTS] ) taken into account and check 
a posteriori that there is consistency between the assumption and the result. 

The thesis is organized as follows: In Chapter 2, we introduce the Proton Spin 
Problem via the polarized deep inelastic scattering experiments and the relation 
between the spin dependent structure function, g\(x), and the flavor singlet axial 
charge. A general review of the chiral bag model is given in Chapter 3 for the next 
discussion. In Chapter 4, we review the axial anomaly and present its contribution 
to the flavor singlet axial charge. Moreover, we show a derivation of the color 
anomaly boundary condition. We address the various static contributions and 
calculate the Casimir effect to the flavor singlet axial charge in Chapter 5. Finally, 
we discuss our result in Chapter 6. 



Chapter 2 

Proton Spin Problem 



Deep inelastic lepton-hadron scattering (DIS) has played an important role in 
understanding the internal structure of hadrons. The discovery of Bjorken scaling 
in the late nineteen sixties provided the basis for the idea that hadrons are made 
up of point-like constituents. The subsequent development of the Parton model 
played an essential role in linking the partons to the quarks via DIS sum rules. DIS 
was essential in the discovery of the missing constituents, identified as gluons, and 
therefore in assembling all different pieces of the hadronic puzzle into a coherent 
dynamical theory of quarks and gluons, Quantum Chromodynamics. 

Polarized DIS, describes the collision of a longitudinally polarized lepton beam 
on a nucleonic target polarized either longitudinally or transversely to an arbi- 
trary direction. It provides a more complete insight into the structure of the 
nucleon than unpolarized DIS. Whereas the latter probes the number density 
of partons with a fraction x of the momentum of the parent nucleon, polarized 
DIS leads to more sophisticated information, namely it determines the number 
density of partons with given x and given spin projection in a nucleon of definite 
polarization. [] 

1 We summarize here the conventions for the Dirac spinors. With four momentum p M = 
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In this chapter, we give a short review of polarized DIS, show the relation 
between the spin dependent structure function, gi(x), and the flavor singlet axial 
current, and discuss some relevant facts about the proton spin problem. 



2.1 Polarized Deep Inelastic Scattering 

In the laboratory frame the differential cross section for the polarized lepton- 
nucleon scattering has the form 

_^--J-^L (2 1) 

dQdE'~ 2Mq* E L » vW ' [ZA) 

where the four momenta of the incoming and the outgoing lepton with mass m 
are k = (E, k) and k! = (£", k'), respectively, and the four momentum for the 
nucleon is P = (M, 0). The momentum transfer is q = k — k' and a is the fine 
(E, p), the Dirac spinors are normalized as; 

v)u = 2E, vh = 2E, uu = 2M, vv = -2M, 

for both the massive and massless case. From these, the following relations can be derived 

u(p)^u(p) = 2p^, v(p)^v(p) = 2p^. 

Incorporating 75 matrix, for a fcrmion of mass M, there is the relation 

u(p, S)n 5 u( P , S) = -v(p, S)n 5 v( P , S) = 2MS», 

with the covariant spin normalized S^S^ = — 1. Additionally, for massless fermion of helicity 
A = ±1/2, the above relation is changed to 

u(p, \)^ l5 u{p, A) = -v(p, X)Xj 5 v(p, A) = lim 2MS»{\) = 4AfA 

M— >0 

Moreover, all states are normalized so that 

(P'\P) = {2irf2E 5 3 (p'-p). 
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structure constant. In eq. (|2.1|) the leptonic tensor is given by 

V = [u(k', s') ltl u(k, s)}* [u(k', s') lu u{k, s)\ , (2.2) 

where s (s') is the spin four vector of the incoming (outgoing) lepton such that 

s ■ k = = s' ■ k' and s ■ s = —1 = s' ■ s' . L^ v can be decomposed into symmetric 
(5*) and antisymmetric (A) parts under fi, v interchange; 

L^(k,s;k's') = L s ^(k, k') + iL* v (k, s; k') 

+L%(k,s;k',s')+tL^(k;k',s'). (2.3) 

Explicitly they become 

L/jj/ik, k) = k^k u + k^ky — g^ v {k ■ k — m ), 
L^(k,s;k') = me^s^k-k'f, 
L%(k,s;k',s') = (k-s')(k^s u + sX-9,uk' ■ s) 

-(k ■ k' - m 2 )(s M s' t/ + s'^Sy - g^s ■ s') 

+(k' ■ s)(s' fl k u + k^s' u - (s • s')(kfj,k' u + k'^k u ), (2.4) 

and 

Lf v = ■me llvaP s' a {k - k'f (2.5) 



where m is the lepton mass |34||. Summation over s' leads to 2L^ U + 2iL flu 



Summation of L^ u over s' and averaging over s gives the unpolarized leptonic 
tensor, 2L^ U . 

Due to the internal structure of hadrons, the hadronic tensor W^ u is unknown 
and is defined in terms of four structure functions as |35j [[3(J; 



W^(q; P, S) = W^(q; P) + iW£(q; P, S) (2.6) 
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with 



— W s 

2M ^ 



—W A 

2M ^ 



P„ 



P-q 



2 %)[ P --^-^ 



W 2 (P-q,q 2 



M 2 



MS^G^P-q.q 2 



+{{P ■ q)S^ - {S ■ q)P^ 



G 2 (P-q,q 2 ) 



M 



(2.7) 



(2.8) 



,s 



(2.9) 



where is the spin four vector of the nucleoli. With these structures eq. ( gjj) 
becomes 

d 2 a 1 a 2 E'r 
dQdE' ~ 2Mlf~E 

The individual terms inside the square brackets can be separately studied 
by considering cross-sections or differences between cross-sections with particular 
initial and final polarizations. Each of these terms is an observable quantity in 
terms of the spin-averaged structure functions Wi, W 2 and of the spin- dependent 
structure functions G\, G 2 . For example, while the unpolarized cross-section 
contains only L^ v W tJ,v,s 



d 2 a unp 
dQdE' 



4 dQdE 1 2M g 4 E ^ 



(2.10) 



the difference of cross-section with opposite target spins contains L A u W^ h 
d 2 a „ „ , , , s d 2 a 



E 



dQdE 



-(k, s, P, S;k , s 
1 a 2 E' 



dQdE 
AL A W^ U ' A 



-(k, s, P, S; k , s 



(2.11) 



2Mg 4 £ ~ ^ 

In the laboratory frame, the cross-section for the inelastic scattering of an 
unpolarized leptons on an unpolarized nucleon, can be written explicitly as 



d 2 a unp 4a 2 E' 2 



dQdE' 



9 9 
2W X sin 2 - + W 2 cos 2 - 



[2.12) 
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where the lepton mass has been neglected. Here 9 is the scattering angle of 
the lepton. This cross section provides information on the unpolarized structure 
functions Wi(P ■ q, q 2 ) and W 2 (P ■ q, q 2 ). In the deep inelastic scattering regime, 
the Bjorken limit is defined by 

Q 2 Q 2 



q 2 = Q' 



oo, 



V 



E-E' 



oo, 



X 



2P ■ q 2Mu 



fixed, 



(2.13) 



lim MWi(P-q,Q< 



and the structure functions obey, so called, scaling for fixed x |37| ; 

= Fi(x), 

lim vW 2 (P-q,Q 2 ) = F 2 (x), (2.14) 

Q 2 ^oo 

Similarly, from eq. ( |2.11|) , eq. ( |2.4|) , and eq. (|2.5| ), the difference of the cross- 
sections with opposite target polarization can be written as; 



dV 



2„s,-S 



d 2 a 



dVtdE' a 


E 

s' 


' d 2 a 


dttdE' 



-(k, s, P, S;k , s) 



8ma 2 E' 
q 4 ~E 



d 2 o 

dndE' (k,s,p, ~ S] k ' }S '\ 

q-S)(q-s) + Q 2 (s-S)^MG l 

)) — 

7 M 



(2-15) 



+Q^(s ■ S)(P ■ q)-(q. S) (P-s)) 

This expression supplies information on the polarized structure functions G\{P ■ 
q,q 2 ) and G 2 (P • q,q 2 )- In the Bjorken limit, they are also known to obey the 
scaling, 

(P ■ q) 2 



lim 



-Gi(P-g,Q' 



lim u(P-q)G 2 (P-q,Q' 



92(x). 



In terms of gi t2 the expression for W^ v can be written as 



W, 



2M 
~P~q 



-[iisafSq 



S ■ q 



SP 9l (x, Q 2 ) + (SP - —± P?g 2 (x, Q 2 ) 

P ■ q 



(2.16) 



(2.17) 
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To get information on the polarized structure functions Gi, G2, we need to 
look at eq. (2.15) with particular spin configurations of the incoming leptons and 



the target nucleons. We consider firstly the case of longitudinally polarized lep- 
tons. The symbol, — > (*—), denotes the spin of the initial lepton along (opposite) 
to the direction of motion and the nucleons at rest are polarized along (S) or 
opposite (-S) to an arbitrary direction S; 

s ^ = -s^ = — (IkLkE), 
m 

with k = , and 
|k| ' 

S" = (0,S). (2.18) 
Choosing the z-axis along the incoming lepton direction, we have 

F = (E, 0, 0, |k|) ~ E(l, 0, 0, 1), 

k'" = {E',k') ~ E'{l,k') 

= E'(l, sin 9 cos0, sin 9 sin <f), cosd), 

S = (sin a cos j3, sin a sin /3, cos a). (2-19) 



This kinematical scenario is depicted in Fig. [2.1| 
Substituting these vectors into eq. ( |2.15| ) yields 

d 2 a^ s d 2 a^- s _ _^E 
dVLdE' dVLdE' ~ Q 2 E 

x \(E cos a + E' cos Q)MG X + 2EE'(cos - cos a)G 2 ] (2.20) 

where 6 is the angle between the outgoing lepton direction, k, and S; 

cos = sin 9 cos sin a cos (3 + sin 9 sin sin a sin /3 + cos 9 cos a 

= sin 9 sin a cos + cos 9 cos a. (2-21) 
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Figure 2.1: The angles defining the kinematical and spin variables of the studied 
polarized cross section are shown 

Here <p is given by <p = (3 — (p. For nucleons polarized along (=>) the initial lepton 
direction of motion or opposite (■<=) to it, that is, a = 0, = 9, eq. ( |2.20|) gives 

d 2 a^ d 2 a^^ Aa 2 E' 



dttdE' 



dVLdE' 



Q 2 E 



{E + E'cos6)MG 1 -Q 2 G 2 



(2.22) 



For transversely polarized nucleons, that is, when the spin of the nucleon is per- 
pendicular to the direction of the incoming lepton, a = 7r/2, and eq. ( |2.20|) yields 

d 2 a^ d 2 a^ Aa 2 E' 2 



sin9 cos <p(MG 1 + 2EG 2 ). 



(2.23) 



dVtdE' dVtdE' Q 2 E 

In the case of ip = tt/2, which corresponds to the nucleon spin being perpendicular 
to both vectors k and k', the difference of cross-sections, eq. (|2.23j ), vanishes. 
The value of such a difference has a maximum when ip = 0, or it, that is, when 
the nucleon spin vector, which is perpendicular to k, lies in the plane determined 
by the two vectors k and k'. 

Experimentally, the polarized structure functions g\ and g 2 are determined by 
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measuring two asymmetries 

_ da^ - da^^ da-** - do-'* 

11 ~ da-^ + da-^ L ~ da^ + da^' { 

where the abbreviation da for d 2 a/dQdE' has been introduced. 

Using the fact that the denominator is simply twice the unpolarized cross- 
section, from eq. ( |2.12| ), eq. (|2.22|) , and eq. ( |2.23| ), the asymmetries become 

Q 2 [(E + E'cos6)MG 1 - Q 2 G 2 ] 



An 



2EE'[2W 1 sin 2 § + W 2 cos 2 § 



Q 2 sm9(MG 1 + 2EG 2 ) 

Aj_ = TTa — 7J-COS09. 2.25 

2E[2W X sin 2 | + W 2 cos 2 § ] y K 1 

It is convenient to write the asymmetries A\\ and A± in terms of the virtual 
Compton scattering asymmetries A\ 2 given by || 

A x = " 1/2 ' ff3/2 , A 2 = 2aTL , (2.26) 

O"l/2+0 r 3/2 0"l/2 + 0"3/2 

where a\/ 2 and 03/2 are the virtual photon absorption cross sections for 7*(1) + 
N{—\) and 7*(1) + N(^) scatterings, respectively, and a TL is the cross section 
for the interference between transverse and longitudinal virtual photon-nucleon 
scatterings. The asymmetries A 12 have the bounds 

l^il < 1, \A 2 \ < y/R, (2.27) 

where R is the ratio of the longitudinal to transverse cross section, R = ai/vr, 
with o T = ((T1/2 + cr3/2)/2. The asymmetries can be written in terms of A\ 2 as 

A\\=D(A 1 +nA i ) 1 A ± = D(A 2 -^A 1 ), (2.28) 

where D is a depolarization factor of the virtual photon, 7] and £ depend only 
on kinematic variables ||34| . The asymmetries Ai j2 in the virtual photon-nucleon 
scattering have relation to the polarized structure functions gi and g 2 ; 

= 9i -I 2 92 = l(9i + 92) 

14 



with 7 = Q/u = Q/(E - E') = 2Mx/^Q z and F 1 in eq. ( gig) . Since in the 
Bjorken limit 7 goes to zero, one obtains 

A,, _ F 2 (x,Q 2 ) A,, 



^2.30) 



D 2x(l + R(x,Q 2 )) D ' 

where F 2 (x,Q 2 ) is the unpolarized structure function in the scaling regime and 
R is given in terms of the unpolarized structure functions 



R 



W2 



1 + 



- 1. 



(2.31) 



Note that the last result of the eq. ( |2.30| ) is from the fact that R can be written 
in terms of the Bjorken scaling functions, 

F 2 (x) 



R 



2xF 1 (x) 



1. 



(2.32) 



in the limit -> 0. 



Experimental results on the polarized structure functions g±(x) for the nucleon 
can be found in the Table f\ 



2.2 The Structure Functions in the Parton Model 

In the parton model the nucleon is regarded as a collection of almost free con- 
stituents, namely the partons, each carrying a fraction x' of the nucleon four 
momentum. Lepton-nucleon DIS can be understood as the incoherent sum of 
scatterings between the lepton and the spin- 1/2 partons M pT7|| . We shall assume 
for our description that the charged partons are quarks and antiquarks, a state- 
ment which was proven historically a posteriori by studying the experimental 
structure function sum rules. The hadronic tensor W^ v can be obtained in terms 
2 We have quoted this table from ref. |3q| . 
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Table 2.1: Experiments on the polarized structure functions gi(x, Q 2 ), g™{x, Q 2 ) 
and gf(x, Q 2 ). 



Exper. 


Year 


Target 


(Q 2 ) 

(GeV 2 ) 


x range 


p target 

= /o5i arslot (^(Q 2 »^ 


Ref. 


E80/E130 


1976/1983 


P 


- 5 


0.1 < x < 0.7 


0.17 ±0.05* 




EMC 


1987 


P 


10.7 


0.01 < x < 0.7 


0.126 ± 0.010 ±0.015 t 


n 


SMC 


1993 


d 


4.6 


0.006 < x < 0.6 


0.023 ±0.020 ±0.015 




SMC 


1994 


P 


10 


0.003 < x < 0.7 


0.136 ±0.011 ±0.011 


© 


SMC 


1995 


d 


10 


0.003 < x < 0.7 


0.034 ±0.009 ±0.006 


@ 


E142 


1993 


n 


2 


0.03 < x < 0.6 


-0.022 ±0.011 


11 


E143 


1994 


P 


3 


0.03 < x < 0.8 


0.127 ±0.004 ±0.010 


© 


E143 


1995 


d 


3 


0.03 < x < 0.8 


0.042 ±0.003 ±0.004 





* Obtained by assuming a Regge behavior A\ oc x 1 for small x. 

t Combined result of E80, E130 and EMC data. The EMC data alone give rf = 0.123±0.013± 
0.019. 
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of the elementary quark tensor w^ u as, 



W(q;P,S) = W^(q;P)+iW^(q;P,S) 

1 f 1 dx' 

= H e l^5 — / T 5 ( x ' ~ x)n q (x',s;S)w fMU (x',q,s), (2.33) 
TT H 2P ■ a Jo X 



q.s 



where n q (x', s; S) is the number density of quarks with charge e q . Here s is the 
spin of the quarks inside a nucleon with the spin S and four momentum P, the 
sum runs over quarks and antiquarks, and x is the Bjorken variable given in 
eq. ( p,13|) . The quark tensor has the same form as the leptonic tensor, eq. ( |2.4|) 
and eq. (|2.5|) , with the replacements fc M — > xP^ and — > xP^ + q^ . After 
summation over the unobserved final quark spin, w^ v becomes 

w^(x, q, s) = w^(x, q) + iw£,(x, q, s) (2.34) 

with the quantities 

w%{x,q) = 2[2x 2 P^P v + xP^q v + xq^P v -x(P-q)g tll/ } 
w^(x,q,s) = -2m g e fl „ a f 3 s a q f3 , (2.35) 

where the quark mass has been taken to be m q = xM for consistency. 

Comparing these equations with the definition of the structure functions 
eq. (|2.8| ), the unpolarized structure functions become 

F i( x ) = \Y. e l < i( x ) 

Z q 

F 2 (x) = x^ej(x) = 2xF 1 (x), (2.36) 
q 

where the unpolarized quark density is defined by 

Q{x) = ^2n q (x,s;S). (2.37) 
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Similarly the polarized structure functions are obtained as 

z q 

g 2 {x) = 0, (2.38) 

where Aq(x, S) is the difference between the number density of quarks with the 
spin parallel (s = S) to the nucleon spin and those with the spin anti-parallel 
(s = -S); 

Aq(x, S) = n q (x, S; S) — n q (x, —S; S). (2.39) 

It is known that in the parton model Aq(x, S) cannot depend on the direction of 
the nucleon spin S, that is, Aq(x, S) = Aq(x) [ j34fl . 

2.3 Relation between the Spin Structure Func- 
tion gi and the Flavor Singlet Axial Charge 

From the previous analysis or from the operator product expansion (OPE) PB| . 
the first moment of the polarized proton structure function defined by, 

T?(Q 2 )= f l £(x,Cf)dx (2.40) 
Jo 

can be connected to the flavor singlet axial current of the quarks by the relation 
T?(Q 2 ) = f gl{x,Q 2 )dx=We 2 q Aq{Q 2 ) = We 2 q {p,S\ q W q \p,S)S", 

JO Z, q A g 

(2.41) 

where Aq is the net helicity of the quark flavor q along the direction of the proton 
spin at momentum transfer — Q 2 . In general, the form of Aq depends on Q 2 . For 
example, in the infinite momentum frame, it becomes 



Aq = / Aq(x)dx = 

Jo Jo 



q\x) + q^(x) — qUx) — qUx) 
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dx. (2.42) 



2.4 The Proton Spin Problem 



At the EMC energies Q 2 < lOJGeV P5|, three light flavors are relevant and the 



first moment of the polarized proton structure function has the form 

T?(Q 2 ) = \ Q Au(Q 2 ) + I Ad(Q 2 ) + 1 A S (Q 2 )) . (2.43) 

In terms of the form factors Q in the forward proton matrix elements of the 
renormalized axial currents ESI, i.e., 



(p, S\ 4 |p, 5) = ia 3 ^, (p, 5| 4 |p, S) = 

(p,S\Al\p,S) = a°S fl , (2.44) 

with 

K = ^7m75 yV>, = ^7^75^, (2.45) 
the sum rule for the first moment is 

T?(Q 2 ) = ^C» s (a s (Q 2 )) (a 3 + ^a 8 ) + (a s (Q 2 ))a°(Q 2 ), (2.46) 

where a s {Q 2 ) is the perturbatively running QCD coupling constant and C\{a s {Q 2 )) 
are first moments of the Wilson coefficients of the singlet (S) and the non-singlet 
(NS) axial currents given by 

2 /n \ 3 



Cf s = l-^i-^(^) - 20.22(^ 



Cf = 1 - — - 1.10( — ) , (2.47) 



1 - 


a s 


43 /a s \ 




Ti 


12 W J 








1 - 




-i.iof^ 




71 


V 7T 



TT 



2 



3 Here we used the normalization of the spin vector of the proton: 



S^S^ = -M 



2 



instead of the previous one. 
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up to 0{a i s ) for three quark flavors. Since there is no anomalous dimension 
associated with the axial- vector currents A 3 and A % , the non-singlet form factors 
do not evolve with Q 2 . The non-singlet form factors are related to the SU(3) 
parameter F and D by 

a 3 = F + D, a 8 = 3F-D. (2.48) 



Their values 50 51 are 



F 

F = 0.463 ± 0.008, D = 0.804 ± 0.008, — = 0.576 ± 0.016 (2.49) 

and from these a 3 = 1.2670 ± 0.0035. From the definitions, the form factors can 
be written in terms of the quark polarizations 

a°(Q 2 ) = Au(Q 2 ) + Ad(Q 2 ) + As(Q 2 ) = £(Q 2 ), a 3 = Au(Q 2 ) - Ad(Q 2 ), 
a 8 = Au(Q 2 ) + Ad(Q 2 ) -2As(Q 2 ). (2.50) 

Before the EMC measurement of the polarized structure functions, a predic- 
tion for rf known as the Ellis- Jaffe sum rule |26|] was based on the assumption 



that the strange sea quark in the proton is unpolarized 

T?(Q 2 ) = ^a 3 + ^a 8 , (2.51) 

without QCD corrections. The measured result of EMC, = 0.126 ± 0.010 ± 
0.015, is smaller than what was expected from the Ellis- Jaffe sum rule: rf = 
0.185 ± 0.003 without QCD corrections and T\ = 0.171 ± 0.006 with the leading- 
order correction. From eq. (|2.48|) , eq. ( [2.50|) , and the EMC result, the quark 
polarizations are obtained as 



Au{Q 2 ) = 0.77 ± 0.06, Ad{Q 2 ) = -0.49 ± 0.06, 

As(Q 2 ) = -0.15 ±0.06, (2.52) 
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and 

AS = 0.14 ±0.17 (2.53) 

at Q 2 = 10.7GeV 2 . The results eq. (gjp and eq. (|233|) reveal two surprising 
things: The strange quark sea has negative non-vanishing polarization, and the 
total contribution of quark helicities to the proton spin is small and consistent 
with zero. These facts raise some puzzles, for example, from where does the 
proton get its spin? why is there negative polarized strange sea quark? how is the 
total quark spin component small? These puzzles are sometimes (inappropriately) 
referred to as the proton spin problem (or crisis). 

The proton spin problem arises from the fact that the experimental results 
seem to be in contradiction with the naive quark-model. The non-relativistic 
SU(6) constituent quark model yields that Am = | and Ad = — |. Therefore, 
from these polarizations one gets AE = 1 and g\ (= a?) = |, which is larger 
than the measured value 1.2670 ± 0.0035 [Q. In a relativistic quark model, the 
quark polarizations Au and Ad are reduced by the same factor of | to 1 and — |, 
and g\ is reduced to | due to the presence of the lower component of the Dirac 
spinor. The reduction of the total quark spin AE to 0.75 requires that the orbital 
angular momentum of the quark, Lq, contributes to the nucleon spin as required 



by the sum rule [52 



l - = IAS + Lq. (2.54) 

Therefore, it is expected that in the relativistic quark model 3/4 of the proton 
spin arises from the quarks and the quark orbital angular momentum accounts 
for the rest of the spin. The MIT bag model, which is a relativistic model with 
QCD confinement incorporated via its boundary conditions, leads to the similar 
value: AE = 2/3. On the other hand, the Skyrme model for the baryons yields 
AE = 0@. 

One way to understand the experimental value AE ~ 0.30, which is smaller 

21 



than the expectation of the quark models, is to introduce a negatively polar- 
ized quark sea. The quark polarization can be decomposed into valence and sea 
components, Aq = Aq v + Aq s . Then, the the total quark spin of the proton 
becomes 

AS = AT, V + AS S = (Au v + Ad v ) + (Au s + Ad s + As,). (2.55) 



The gluons can induce a quark sea polarization through the U(1)a anomaly |28 



which cancels the spin from the valence quarks when the gluon has negative spin 



component f38] . 



Another way is to use the axial anomaly directly in calculating the flavor 
singlet axial current. In other words, the experimentally measured quantity is 
not merely the quark spin polarization AS but rather the singlet form factor 
(the flavor singlet axial charge), to which the gluons contribute through the axial 
anomaly as 

a°(Q 2 ) = AS - N F ^r^Ag(Q 2 ), (2.56) 

Z7T 



where Ag is the polarization of the gluons and Np the number of flavors p4 
These explanations, and possibly others, could be reconciled if one were to es- 
tablish that they are gauge dependent statements, while the measured quantity 
is gauge- invariant p>l. Incorporating the gluons, the spin sum rule becomes 



i = l - AS + L Q + Ag + Lq (2.57) 

with the orbital angular momentum of the gluon, Lq, and the integral of the 
polarized gluon distribution, Ag. 

The analysis of the flavor singlet axial charge and the gluon spin in the chiral 
bag model will be discussed and compared with those of the MIT bag model in 
Chapter 5 after introducing anomalies in Chapter 4 . 
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Chapter 3 

The chiral bag model 



In this chapter we give a general overview of the chiral bag model as initially pre- 
sented [Q. We review its definition in terms of quark, gluon and meson degrees 
of freedom. Here we shall be dealing with a Lagrangian that is classically gauge- 
invariant. We discuss the solution with this Lagrangian obtained by using the 
hedgehog ansatz, including the effects on the vacuum structure. Our discussion 
will incorporate the rf meson since it will be relevant for later purposes. It turns 
out that due to color anomaly, this theory is not gauge invariant at the quantum 
level. We avoid in here the complications arising from the quantum structure of 
the theory, relegating this subject to the next chapter. 



3.1 Model Lagrangian 

The chiral bag model is a field theoretic description of hadron structure whose 
aim is to represent QCD in the low energy regime. This description separates 
space-time into two regions by a surface, the bag, in which different effective 
realizations of the underlying theory, QCD, are used to represent the dynamics. 
The bag, which is closed in space, defines an interior region, conventionally called 
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quark phase, which is described by means of quark and gluon fields. The exterior 
region, called mesonic phase, is defined by an effective mesonic field theory in 
accord with the requisites of Weinberg's unproven theorem |T7|]. The bag, the 
surface separating the two phases, serves to connect the two types of degrees of 
freedom through boundary conditions, whose structure resembles the bosoniza- 
tion relations in two dimensions. 

The motivation for this sophisticated description lies in the properties of the 
fundamental theory, QCD, which the model implements in a dynamical fashion. 
Let us be more precise: 

(1) Color Confinement: the bag is responsible for confining the color degrees 
of freedom (quarks and gluons) and the boundary conditions on it implement the 
non perturbative character of this property. Despite the apparent weak interac- 
tion between these fields in the quark phase, the fact that they are represented 
by cavity modes satisfying the boundary conditions, confers them a non pertur- 
bative character very different from that of free Fock states, even in the case of 
an empty mesonic sector |§. 

(2) Asymptotic freedom: It is known that quarks and gluons in QCD interact 
very weakly at large momentum transfers, i.e., short distances. This important 
property of the theory, associated with the negative sign of its (3 function |57] 
is responsible for the slow logarithmic deviations from scaling in deep inelastic 
scattering. In the bag description it is implemented by the perturbative treatment 
of the interaction between quark and gluons in the interior region. 

(3) Spontaneous broken chiral symmetry: Nature, and therefore QCD, real- 
izes chiral symmetry in a spontaneously broken fashion, i.e. the flavor symmetry 
SU^Nf) x SUft(Nf) of the currents is broken down to SUy(Nf) by the vacuum. 
This phenomena is implemented by the mesonic effective theory outside, which 
incorporates the pseudoscalar mesons, the required Goldstone bosons in the chiral 
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limit. Through the boundary conditions this phenomenon transfers to the inte- 
rior. The ultimate objective of the model is to encompass both long-wavelength 
and short-wavelength regimes, with the Cheshire Cat principle defined below 
bridging the two regimes. 

In this well defined scenario with a given chiral bag Lagrangian, the boundary 
plays a crucial role because it relates the degrees of freedom of the two phases in 
a manner which preserves all the symmetries and their realization. 

The chiral bag model as described above can be implemented by the following 
Lagrangian density 



where Cq, Cm and Cqm describe the dynamics for the quark and gluon fields 
inside the bag, the meson fields outside, and the interaction between the quark 
and meson phases at the bag surface, respectively. Here O^, which is needed 
to define the quark phase inside the bag only, it gives 1 inside the bag and zero 
outside, Q B = 1 — @b> and the bag delta function A B is defined by A B = 
—n^d^Q B where represents the outward normal unit four vector. B is the so 
called bag constant and corresponds to the energy density required for creating 
the bag in the QCD vacuum. 

The Lagrangian density for the quark phase in case of SU (3) flavor symmetry 
is given by 



£ = (£q — B)Q B + £>m@b + £qm& b , 



(3.1) 



£q = ${i-f\(% -X)~M)i>- g s G a ^Y^ ~ \ G %G a ^ (3.2) 



with 



G 




(3.3) 
(3.4) 
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Here ip represents the quark fields, the gluon fields, M the current quark mass 
matrix m u ~ ~ 0, m s m 150 MeV. Asymptotic freedom is realized by allowing 
the interaction between the quark and gluon fields to be treated perturbatively 
with respect to the effective bagged QCD coupling constant g s . The Lagrangian 
has the same form as that of QCD, but it is only meaningful in the weak coupling 
regime. The quark field is arranged into the fundamental representation of flavor 
SU{3) 

V>= d . (3.5) 

V s J 

A" (a = 1, 2, • • • , 8) are the Gell-Mann matrices associated with color and we use 
the normalization TrA^A* = 28 ab . 

The meson phase is described by the following Lagrangian density 



-aTr(M(U + - 2)) Ir^ml, ( Tr(ln[/ - lntf)) , (3.6) 

16N F 1 V / 



where the chiral field U is of the form 

./() In 



U = exp ( i^r + i^X) (3.7) 



with /o = J^U, and 



2a = (uu) = (dd) ~ (ss) Q 

m lfl _ ™ 2 K fl 



(3.8) 



m u + m d m s 

where A a (a = 1, 2, ■ • • , 8) are the Gell-Mann matrices associated with the flavor 
symmetry in this case, f n is the pion decay constant^], m n , ihr- and m v ' repre- 
sent the pion kaon and rf masses, respectively. The rf has been introduced by 

Unless otherwise specified we assume /„• ~ fx- The difference appears at higher order in 
the chiral counting. 
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extending the chiral field U from SU(3), as appears in the chiral hyperbag |JT 



to U(3). In this way the r( is decoupled from the other pseudoscalar mesons. It 
plays an important role in the flavor singlet axial charge []58] due to its flavor 
singlet structure. The second term in the above Cm is the one introduced by 
Skyrme to stabilize the embedded SU(2) soliton solution. From an analysis of 
the nucleon axial form factor qa in the Skyrmion model |12|] , the parameter e can 
be fixed to 4.75 throughout whole bag radius. 

The last term, the so called Wess-Zumino-Witten term Cwzw @, comes 
from the requirement that an effective theory should have the same symmetries 
and anomalies as the fundamental theory, at its validity scale. Its explicit form, 
which can only be written as an action, is 
N r 

Twzw = -i^h / d " x ^ Xpa MU ] d^d v UU^d x UU^d p UU^d a U), (3.9) 
24uir z J m 

where the integral is defined on the five dimensional manifold M = B x S 1 x [0, 1] 
with B, the three-space volume outside the bag, and S 1 the compactified time. 
The extension [0, 1] is needed to be able to write the Wess-Zumino-Witten term 
in a local form. 

The interaction between the quark field and meson field on the bag surface is 
given by 

£qm = ~^U 5 iP = --(iPlUiPr + ^jrUHl) (3.10) 



with tjj^L = |(1 ± 75 ) and U 5 = exp m^T/Afo + A • 7r// 7r )J. This interaction 
provides quark confinement classically in a chirally invariant way. Note that no 
interactions between the gluon field and meson field appear at the classical level. 

In case of massless quarks, this Lagrangian is invariant under flavor SUl{3) x 
SUr(3) transformations. Noether's theorem gives the following conserved cur- 
rents 



j;> r = ~^(i+7 5 )7^©b 
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+ 



4 lbe z 



N / \ a 

+ 48^ wTr IT U ] d u Utfd p UlP& r U 



Qb, (3-11) 



J? L = ^(i-7 5 )7^e 



+ 



+ 



~Tr{\ a d^U) - -^{[\ a ,cTUU^[d,UU\d v UU^]) 
N c „ /A a 



Tr [ —d u UU*d p UU*d a UU^ ] 



487r 2^P--^y- - - - - yj ^ ( 3 - 12 ) 

where the index a runs over 1, • ■ ■ , 8. The vector and axial vector currents can 
be constructed from these currents as 

v; = r, R +J a , L , 

A * = j;> R -r+. (3.13) 

The baryon number current corresponding to the t/y(l) symmetry of the La- 
grangian is 

= in^Q B + ^e^ Xp Ti(U^d u UU^d x UU^d p U)& B . (3.14) 

The last term corresponds to the topological winding number arising from the 
Wess-Zumino-Witten term after proper gauging \\59\ . The conservation of this 
term is a consequence of topology in case of the SU(2) symmetry. The Ua(1) 
symmetry of the Lagrangian yields a flavor singlet axial vector current of the 
form 

4 0) = ^757^©b + 2f n d,r ] '(x)Q B , (3.15) 
which is broken through the well known axial anomaly |28| providing the rj' with 
its mass. Its role in the proton spin problem will be discussed in chapter 5. 

The Hamiltonian, of the chiral bag model, can be derived from the Lagrangian 
and turns out to be 

H = [ al 3 r ^( - ia • V + p)ip + - I d 3 r ipU 5 ip 

JB V / 2 JdB 

28 



r r f 2 
+ d 3 r ^ 
is L 4 



-TrCSblT+SbZ/ + <9^[/) 



+crTr(UM + M ] U ] - M - M ) + H wzw . 



(3.16) 



Finally, we should mention the magnitude of the quark-gluon coupling con- 
stant g s , a parameter of the model. Although g s has a scale dependence, governed 
by the f3 function, it is customary to take it at a fixed value throughout the vol- 
ume, as corresponds to lowest order perturbation theory. The MIT group fixed 
g s by studying the masses of the proton and the delta in their model 0, and 
obtained a s = g 2 /47r = 2.2, which we will use for computational purposes. 



Our aim here is to describe a solution to the equations of motion of the model 
Lagrangian. Our first approximation will be to consider the bag as a static sphere 
of radius R. The symmetries of the Lagrangian are instrumental in finding the 
adequate solution. In the quark sector we recall that the up and down quark 
masses can be neglected since they are small, the strange quark mass is neither 
so small to be neglected nor so large to allow a heavy quark treatment. Moreover, 
the strange quark mass breaks not only chiral symmetry SUl(3) X SUr(3) down 
to SU V (3) but also the symmetry SU(3) down to SU V {2) x U Y (1), i.e., SU L (3) x 
5^(3) - SU V (2) x U Y (1). 

The symmetry breaking scheme of the quark phase suggests that the meson 
phase can be described by the classical configuration of the chiral field Uq which 
is the SU{2) hedgehog solution embedded in SU{3). The explicit form of Uq is 



3.2 The Hedgehog Solution 
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given by [60 

/ e tT-M(r) o 

U = exp{iX i f i B(r)) = q ^ ), (3.17) 

where r» are the Pauli matrices. 9(r) is called the chiral angle. We do not include 
now the rf meson in our description, but we will do so later when its presence 
becomes relevant. 

From the Lagrangian of eq. Q3.1| ) and the above hedgehog ansatz, if the quark- 
gluon coupling is turned off, the quarks, in the spherical cavity approximation, 
satisfy the following equation of motion and boundary condition 2 

ij^iP\ r=R = 0, r<R, (3.18) 
in M YV|r=Ji = Ufy\ r=R . r = R, (3.19) 



The ansatz eq. ( |3.17| ), the equation of motion eq. ( |3.18| ), and the boundary con- 



dition eq. ( |3.19| ) for the quarks show that the strange quark is decoupled from 
the u and d quarks. In addition, due to the Pauli matrices of SU{2) flavor space 
in the hedgehog ansatz, the u and d quarks form the multiplet of the grand spin 
operator K defined by 

K = J + I, (3.20) 

where J is the total spin and I is the isospin. This is called the hedgehog quark 
state. Denoting the wave functions of the hedgehog quark state and the strange 
quark, respectively, by <^(r)e~* £n * and (j)^(r)e~ luJnt with the appropriate quantum 
number n, they satisfy the following equations of motion and boundary conditions 

-za-V^(r) = e„y£(r), r < R, 

-i7-?^(r) = e^^(r), r = R, (3.21) 



2t 



Eq. (3.19) yields the quark confinement condition, iip'y ■ rip = 0, as can be seen by using 



that eq. ( |3.19] ) is changed to iip~f ■ f = ipUs under the hermitian conjugation. 
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-ia-V^(r) = u n (p a n (r), r < R, 

-n-r<p s n (r) = <(r), r = i2, (3.22) 

where 9 S = #(i?) and 
have been used. 

Incorporating the color degrees of freedom \a), the quark field ip can be ex- 
panded in terms of these wave function as 

1>(r,t) = £ <P h n (r)e-^\a)a a n + £ ^(r)e fa »«|a># 

71,0;, £"71 >o <o 

+ E C(r)e-^|a><+ £ C(r)e^V)< f - (3-24) 

n,a, w„>0 n,a, w n <0 

a (6+) is the annihilation operator for the positive (negative) energy hedgehog 
quark and c (eft) the annihilation operator for the positive (negative) energy 
strange quark. The operators, a, b, c, d, satisfy the usual anti-commutation rules: 

R f , O = {K\ €) = K\ ci) = K*, c'} = Smn6 aa > (3.25) 

vanishing all other anti-commutators. The quark vacuum is defined by a"|0) = 
6"|0) = c"\0) = d"\0) = 0, that is, all the negative energy eigenstates of the 
hedgehog and the strange quarks are filled with three different colors. 

The hedgehog quark state has the following quantum number: K, the grand 
spin such that the eigenvalue of K 2 is K(K + 1), M K) the eigenvalue of third 
component of K, P, the parity, and finally n, the radial quantum number. It is 
convenient to introduce additional quantities such as k = P(— 1) K and e, the sign 
of energy eigenvalue which makes the radial quantum number a positive integer. 
The hedgehog quark state will be denoted by |m), i.e., m denotes the set of 
indices {K, M K , P, n, k, e}. 
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The eigenstate \K,Mk) of K 2 and K z can be constructed by a linear com- 
bination of the eigenstates of the total spin operator and the eigenstates of the 
isospin operator. With the help of the eigenstates of the total spin operator 
J = L + S, there are four combinations for K^O. Because of the parity, in terms 
of \K,Mjc)i given in the appendix, the wave function for the hedgehog quark 
state can be written as 



'i) for k = +1 



Vm = a Ni . a . 

%a- ■ rj K+1 (e n r) 

+ l\ 2 [ . 3K } Snr \ \\K.M K h. (3.20) 
-icr ■ rj K _i(e n r 




and 



(ii) for k = — 1 

Vm = aN 1 \ _ a I \K,M K ) 3 

-icr ■ rj K (e n r) 

+ l\ 2 \ . jK : l{£ f\ \\f<^):- (3.27) 
icr ■ rj K (e n r) 

where jK{x)'s are the spherical Bessel functions, N\, N% the normalization con- 
stants: 



( 1 ^r. 



\W> Q n (f K (Q n ) + f K+1 (Q n )) - 2(K+l)j K (Q n )j K+1 (Q r 



at _ (J: ^ \ 2 co 

with f2 n = e n -R, and a, (3 constants with the condition, a 2 + f3 2 = 1. 

Substituting these wave functions into the boundary condition eq. fl3.21|) , the 
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energy eigenvalue e n and the constants a, (3 are determined by the linear equation 
' «(l-H^)^+i-cos^ -sin^Hvggp^ > 
^ - sin s 2 ^2K+t 1] 3k+i k(i + j0^jjK-i + cos 9 s j K J 

= 0. (3.29) 

The energy eigenvalues are obtained from the determinant of this matrix which 
is of the form 

COs6» s {j|(fi n ) - jK+l(tin)jK-l(ttn)} - Kj K (tt n ){ j K +l(^n) ~ Jjf-l(^n) } 

|Sin S'|(fin) = 0, (3.30) 




for an arbitrary K > 0. The equation for if = is obtained simply by setting 
Jk-1 = 0: 

Kcos6>,Ji(fi n ) - (1 + Ksm9 s )j (Q n ) = 0. (3.31) 



In Fig. |3.1| the lowest energy level is drawn as a function of the chiral angle 9 S 
for K = of positive parity. By the structure of the eq. (|3.30|) , the energy levels 
are degenerate with respect to the quantum number M#-. For arbitrary 9 S , the 
energy spectrum is asymmetric with respect to e n = 0. The energy spectrum 
becomes though symmetric for specific values of the chiral angle, 9 S = nn and 
9 S = nir + 7r/2 with (n = 0, ±1, ±2, • ■ •). The case 9 S = nir corresponds to 
the MIT bag model and all quark states have partners for negative energy. For 
9 S = mr + TT/2, the energy spectrum is symmetric except for the zero-energy state. 
Furthermore, the symmetry of the Dirac equation and of the boundary condition 
give the energy spectrum the following symmetries 

e m (e 8 ) = e n (ir + e 8 ), (3.32) 
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Figure 3.1: The lowest energy level of the hedgehog quark state for K = 0. 



where m = {K, Mk, P, n, K, e} and n = {K, Mk, —P, n, —K, e}, and 



(3.33) 



where m = {K, Mk, P, n, k, e] and n = {K, Mk, —P, n, —k, — e}. 

From the Lagrangian for the meson phase with the hedgehog ansatz (classical 
configuration), eq. (|3.17|) , 



U (r) = exp(i\i ■ fi9(r)), 



(3.34) 



the equation of motion for 6{r) and the boundary condition are 

'd6\ 2 2sin 2 



2sin 2 fl\d 2 # 2d6 sin2# 
1 + s— -r^ + 



r 2 / dr 2 r dr 



(r + 2 sin o) — 

e dr 



V 

dr J 
for r > R, 



-~(H| / d 3 r ip^f -T^-hj) |H) , 

2 V&B 



for r = R. 



(3.35) 



(3.36) 
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Substituting the hedgehog ansatz into the Hamiltonian eq. ( |3.16| ), the contribu- 
tion from the meson phase becomes 



2tcU r 2 A { (d9\ 2 sin 2 9 sin 2 9 
E mes0 n = / r dr { ( — ) + 2 — — + 



Ft. \\drj r 2 r 2 



d9\ 2 sin 2 9 



+ 



(3.37) 



where the mass terms have been omitted. Since for the static case Em = —Lm, 
the minimization of the energy with respect to the variation of 9(r) leads also to 
the equations of motion. 

The strange quark has the same energy spectrum as that of the MIT bag 
model since there is no classical configuration associated with the kaons. The 
strange quark states |m) are described by four quantum numbers, j, the total 
spin such that the eigenvalue of J 2 is j(j + 1), rrij, the eigenvalue of the third 
component of the total spin, P, the parity, and n, the radial quantum number 
and, for convenience, the two indices, k, which has the value ±1 corresponding to 
j = I ±1/2, and e, the sign of the energy eigenvalue, i.e., m = {j, rrij, P, n, k, e}. 

Using the eigenstates \j,rrij) K of the total spin J appearing in the appendix 
as basis, the wave functions for the strange quark become 



(i) for k = +1 




icr ■ rj l+1 (u n r) 



\j,mj) K=+1 , (3.38) 



and 

(ii) for k = — 1 




-ia ■ rj l _i(a; n r) 



|j,m j ) K= _ 1 , (3.39) 



where Ni and N2 have the same form as those for the hedgehog states replacing 
Q n and K by X n (= u n R) and I, respectively. Using the boundary condition 
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eq. ( |3.22| ), the energy eigenvalues are obtained from 

jl+1 (X n ) = + 3l (X n ), (3.40) 

for k = +1, and 

ji-i(X n ) = -ji(X n ), (3.41) 

for k = — 1. The energy spectrum is degenerate with respect to the quantum 
number rrij and has the property 

X m = ~Xn (3.42) 



due to the invariance of the Dirac equation and the boundary condition eq. ( 3.22|) 
under the CP-operation. Here m = {j, rrij, P, n, k, e} and n = {j, rrij, —P, n, —k, — e}. 

3.3 The Baryon Number Fractionization 

As can be seen in Fig. |3.1| , the energy level dives into the Dirac sea at 9 S = — 7r/2. 
This means, that even if one has unit baryon number at 9 S = by putting 
N c quarks into the bag, there can be a leakage of baryon number as 9 S varies. 
This fact brings to baryon number fractionization when only the quark phase is 
considered. 

By filling up the valence quark states and all negative energy eigenstates 
with quarks of different colors, the K = ground state hedgehog baryon can be 
constructed as 



cki ■ ■ -a 



a'lO), for -| <fl s <0 



|H )H ^ ° QJVc '~" 2 -" s " (3.43) 

|0), for - 7T < 6 S < -f , 

where |0) is the vacuum state defined by a"|0) = 6"|0) = c"|0) = d"\0) = and 
the subscript "0" in the creation operator indicates that K — 0. The baryon 
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number is given by 

B B = ^-(E \ [ d^Vty |H ). (3.44) 

JM C JB 

Applying Wick's contraction leads to 

f 1 + B vac for - ? < 9 S < 0, , 
B# = < 2 - ' 3.45 

1 B vac for -7r<^<-|. 

While the usual vacuum cannot carry any baryon number, the quark vacuum 
of the hedgehog state can get an induced baryon number through a non-trivial 
polarization by the interaction with the meson phase outside the bag as first 
pointed out by Vento et al [|6T|. The induced baryon number of the hedgehog 



quark vacuum can be obtained by evaluating the regularized spectral asymmetry 



B vac = Km ( - ^ E sign(^)e rM ) 



±(0 s -sin0 s cos0 s ) -|<0,<O, 
l + i(0 s -sin0 s cos0 s ) -7r<0,<-f, 



(3.46) 



where r in the first line is introduced for regularization. The non-trivial vacuum 
polarization and the non- vanishing baryon number of the hedgehog quark vacuum 
result from the CP-symmetry breaking in the energy spectrum. 

As discussed in the previous section, the baryon current, eq. (|3.14j) , gets a con- 
tribution from the meson phase outside the bag through the Wess-Zumino-Witten 
term 0. Substituting the classical configuration (hedgehog solution), eq. ( |3.17| ), 



into the meson part in eq. ( 3.14 ) yields the baryon number 



B 5 = J_d 3 re 0ljk Ti(U^d l UU^d 1 UU^d k U) 



— (6 S -sin0 s cos0 s ), (3.47) 

7T 



3 That is for Nf = 3. For Nf < 3, the argument is indirect as discussed in the previous 
section. 
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where the condition 9(r — > oo) = has been used. Therefore, to get baryon 
number one, two quantities, B B and Bg, should be added. In other words, 
although the quark is confined classically, the quantum fluctuations due to the 
hedgehog solution of the meson induce a leakage at the bag surface of baryon 
number so that there is a contribution to the baryon number from the outside 
region. That is 

B = B B + Bs = l. (3.48) 

This mechanism is known as the baryon number fractionization, and, as just seen, 
in the SU (3) case it is identical to that previously studied for the SU (2) model 
& 



3.4 The Cheshire Cat Principle (CCP) 

We have seen that the baryon number in the chiral bag model (CBM) arises from 
summing the contributions arising from the quark and meson phases. This is a 
particular case of a general statement, namely that for any observable O in the 
CBM its value arises from adding the contribution of both phases, 



O = O b + O b , (3.49) 

where Ob and Ob are the contributions from the quark and meson phases, re- 
spectively. It may be recalled, that in the case of the baryon number, although 
each contribution independently depended on the bag radius R, its sum did not. 
This is because baryon charge is a topological quantity. It has been observed 
that this is a general trend, i.e., when a correct calculation for any observable 
is performed within the CBM, the result tends to be almost radius-independent 
over a sizeable range of R |63| . This statement has become known as the 
approximate Cheshire Cat Principle (CCP) |L5 |. 
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In order to understand the profound meaning of the CCP we have to recall 
some results from quantum field theories in 1+1 dimensions. In this case fermionic 
theories are exactly bosonizable, i.e., one can write for any fermionic theory a 
bosonic theory which leads exactly to the same S-matrix. Thus in 1 + 1 dimensions 



the Cheshire Cat Principle is an exact statement and its meaning very clear ||64 |. 
Let us divide space into two arbitrary regions. In one of them we describe the 
physics by means of a certain theory of fermions. In the other by its equivalent 
bosonic theory. The boundary conditions, which couple the two theories, arise 
from the bosonization rules associated with given symmetries. Any observable 
one calculates arises from the addition of the contribution of the two sectors and 
naturally it is independent on the position of the boundary. Thus the Cheshire 
Cat Principle is a corollary of exact bosonization and the proper definition of the 
boundary conditions. One can phrase this freedom in terms of a gauge symmetry 
In four dimensions there is no exact bosonization technique known up to 
date. This is because one would in principle need infinitely many mesonic degrees 
of freedom to write a theory equivalent to a fermionic theory. Thus the CCP 
can, in general, be only an approximate statement. The exact CCP for the 
baryon number is a special case because of its topological character, i.e., from 
all possible mesons fields only the hedgehog carries baryon number. Therefore 
the CCP transforms from a corollary of exact bosonization in 1+1 dimensions 
to a predictive statement in 3+1 dimensions. It basically asserts the quality and 
indicates the limitation of our effective theories and calculations. The closer our 
theories represent the true theory in their corresponding regime and the better 
we perform our calculations, the larger will be the range of radius independence 
of our observables. 

Many calculations have been performed for different observables and in all of 
them a certain degree of radius independence has been observed [J65|. We will 
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show in this thesis the realization of the CCP in a very complex physical scenario. 



3.5 Collective Coordinate Quantization 

We have studied the hedgehog solution for the ground state of the baryon in the 
previous section. This solution does not carry spin nor isospin and therefore does 
not correspond to any baryon of the spectrum. It can be regarded as a super- 
position of physical B = 1 baryons with various spins and isospins constrained 
by the relation K = J + I = [103]. In other words, the independent spin 
and isospin symmetries of the baryons are mixed up in the K-symmetry of the 
hedgehog solution. By using the collective coordinate quantization method, this 
problem can be overcome and baryons with the appropriate quantum numbers 
can be obtained. 

We shall consider firstly the SU (2) case and then its extension to SU(3) which 
contains the Wess-Zumino-Wess term . 

The hedgehog solution is degenerate in energy with respect to an arbitrary 
constant rotation in 577(2) space A, which transforms the fields as 

ijj _> ijj' = Aip, 

U -> U' = UAU\ (3.50) 

and which can be parameterized as 

A = a + ia. ■ r, (3.51) 

where the parameters are constrained by a 2 , + a 2 = 1. Allowing A to be time- 
dependent introduces three independent collective coordinates. Substituting these 
new fields into the Lagrangian leads to 

L' = L Q + LZ™ k + LZ son i (3.52) 
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where Lq is the original Lagrangian and 

I Jb 

^meson _ —X mcson UJ 2 . (3.53) 

Here u) represents a rotational velocity, which is defined by 

u> = — ztr [tA<9 A] = a a — ad + a x a (3.54) 

and X meson a moment of inertia arising from the meson phase due to the collective 
rotation, which is given by 

'3tt 



roo 

J,,,,..,,,, { — r 2 dr sin 2 9 



/2 | l/dey | sm^ 

71 e 2 \dr J r 2 



(3.55) 



Note in L™f on , eq. ( |3.53D , that the rotational effects associated with the hedgehog 



appear in the mesonic sector to second order in u>, while in the quark phase in 
first order in u). The equation motion for the quark changes to 

^-yT-wj^ = (3.56) 

with the boundary condition given by eq. (|3.19|) . With the adiabatic assumption, 
i.e. slow rotation, the additional terms in the equation of motion can be treated 
perturbatively. The single quark eigenstate obtained by means of standard time 
independent perturbation theory is given by 



Mrl = ^ + \ E (m l T l1 n W ) + " " " , (3.57) 

m ^ 777, r 



where ip$ represents the n-th eigenstate of the unperturbed equation with eigenen- 
ergy e^ m . The baryon state is also modified by the well known Thouless formula 

|H)=exp( £ \- {p [l^ h) ■ ala h \R )\ (3.58) 
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where |H ) is the unperturbed baryon state given by eq. ( |3.43| ) and a) p (ah) is 



the creation (annihilation) operator for a particle (hole) state. The energy of 
the system arises from both phases. The contribution from the quark phase can 
be calculated by taking the expectation value of the Hamiltonian operator with 
respect to the baryon state eq. ( |3.58|) ; 



(H| J B d 3 r ^ol ■ iV^ |H) 



(H|H) 
1, 
2' 



-^quark + 7^quark<^ 2 H , (3.59) 



where the moment of inertia from the quark phase, X quar k, is defined by 

1 \(p\r z \h)\ 2 

Xquark = - }^ Q _ n ■ (3.60) 

p£H ,/ieH fc P t h 

Here, only r z appears because we are choosing the axis of rotation along the in- 
direction. The contribution from the meson phase is obtained by substituting U' 
into the meson part of the Hamiltonian eq. ( |3.16| ); 

-^meson -^mcson 2"^ meson ' t ' ' (3.61) 

where -E^ eson and X meson are given in eq. (|3.37|) and eq. Q3.55|) , respectively. In- 
cluding the volume energy, the energy of the baryon becomes up to second order 
in oj 



-^baryon — ( -Xquark + -^mcson + g 71 ""^ ^) ~^ 2 (~^ c l uark ~^ -^nieson) w 
1 



= E + -lu 2 . (3.62) 

Let us proceed to describe the Isospin and Spin in this formalism and see how 
they enter into the energy expression. Substituting the fields given in eq. (|3.50|) 
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into the expressions, eq. ( |3.12| ) and eq. ( |3.13| ), of the vector current, replacing A a 
by r in the SU(2) case, and integrating the time component of the current over 
space, we obtain the isospin in terms of the collective variables 

Tj = g ( X ^ r ^ Ti ^ ) + ImcsonUJi 

= (^quark "I" ^meson ) = l^ii (3.63) 

where the expectation value for the quark fields has been taken with respect to 
|H). 

The conjugate momenta I1 M (/i = 0, 1, 2, 3) associated with the collective co- 
ordinates (/i = 0, 1,2,3) can be derived from the Lagrangian, eq. ( |3.53|) , and 
yield 

9L 

n = — = 2T-a, 

n = 5 = 2(-Td + axT). (3.64) 
da 

These relations lead to an expression for the isospin in the form 

T = - f alio - a n + a x IlV (3.65) 



2 

Requiring the commutation relation 

[<V,n„] = iSfu,, (3.66) 

the quantum mechanical isospin operator can be represented as 

i ( d d d \ 

Ti= A ao ^r ai d^~ 6ijkaj d^)- (3 - 67) 

We now proceed with the description of spin. For K = the space rotation 

turns opposite to the isospace rotation, therefore, the quantum mechanical spin 

operator can be built by replacing a with its negative value 

i ( d d d \ 

Ji= 2l Oi ^- a0 ^- 6w *^J- (3 ' 68) 
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With these results and the fact, [J, T] = 0, the energy of the baryon can be 
written, keeping in mind that K = 0, as 

T 2 J 2 

-^baryon = Eq + — = E + — , (3.69) 

This expression may be interpreted as the energy of the rotating top with the 
moment of inertia X. Since we are concerned with the static case, this energy can 
be regarded as the mass of the baryon. In terms of the eigenvalues of the isospin 
and the spin, the mass becomes 

M baryon = M + = Mo + J(J 2 ^ 1} , (3.70) 

so that the corresponding masses of the nucleons and the As are 

M N (I=J= 1/2) = M + ~, 

M A (I = J = 3/2) = M + ^. (3.71) 

The extension to SU(3) case is more complex. The hedgehog is sitting in 
SU(2) and only the collective coordinates are extended. By replacing A in 
eq. ( p.50D by a SU(3) matrix, eight collective variables q a (a = l,---,8) are 
defined through the relation 

l -\ a q a = —A^d A. (3.72) 

Substituting the fields Aip and AUA^ into the Lagrangian eq. (|3 . 1| ) , it becomes 
in terms of the collective variables 

L = L Q {A, A) + L M {A, A) + L QM - ^ttR 3 B, (3.73) 



with 



Lm = -ZmesonQi + ~Z meson Qm + 2y/3^^ 8 " ^™ eson, 



Lqm = -77 / d 6 r m^B- (3.74) 

Z JOB 
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Here the static meson field U(r) = exp(i\ifi8(r)) has been integrated out to 
get the L M . The index i (M) in the L M runs 1, 2, 3 (4, • • • , 7). The quantities 



-^meson j 1" meson i B 'g , and iWj 





meson 



are given by 



In 
T 



NT "Vdr sin 2 #|/ 2 + i 
r I e 2 

/•oo 

2tt / r 2 dr(l 

JR 



d6\ 2 
drj 



sin 2 9 



COS I 



^ + 4e 2 



dfl\ 2 2sin 2 fl 
(ir/ r 2 



(3.75) 



B, 



meson 



The X n 



[0 S — sin 6 S cos S ], 

7T 



roo [ f2 

4./ r A i 



d#\ 2 2sin 2 fl 
dr/ r 2 



+ 



sin 2 fl 

2e 2 r 2 



fd6\ 2 
\dr) + 



sin 2 fl 



and X' r 



are the meson contributions to the moments of inertia. 



Note that in case of R = 0, they are just the moments of inertia of the standard 
SU(3) Skyrmion |12] [68|. The Bg is the baryon number carried by the meson 
arising from the Wess-Zumino-Witten term. Since the change of the meson field 
by the collective rotation occurs to second order in the time derivative of A, 
9{r) has been substituted by 6o(r), the value which minimizes the energy(mass) 
M meson ineq. (&7§. 

For the quark field, the effect of the collective rotation appears in first order, 
as did in the SU (2) case. Replacing t-w in the SU (2) case by \ a q a gives the 
equation of motion in the form of 



and the boundary condition 



i if ■ yA,//,, )r = 



(3.76) 



(3.77) 



Assuming an adiabatic collective rotation, the change in the single quark eigen- 
state can be calculated by standard time-independent perturbation theory. Tak- 
ing the wave functions of eq. ( |3.23|) as the unperturbed solutions, the single quark 
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eigenstates are given by 



m = ^ +E «z« (r1+E vi^ +0( . 2)i 
<m = + e N ~ _ ,o > ( 3 - 78 ) 

where m denotes the hedgehog state and n the strange state, respectively. As 
mentioned, i = 1,2,3 and M = 4,5,6,7. Note that Ag does not contribute to 
these perturbations of the wave functions because 0^ and 0^ are the eigenstates 
of A 8 with eigenvalues 4= and — 4=, respectively. The effect of A 8 is only to shift 
the eigenvalues of the quark fields. The perturbation modifies the ground state 
in a form analogous to eq. ( |3.58|) 

im - exD^ V 1 (zalAigila) t 6 t , v 1 (m\X M g M \n) 

|n; - exp i ^ 9 _ c-o a ™°™ + 2^ 9 P o _ , ,o a ™ a n 

x - 1 (m|A M gAf|2l) t ht _i_ V* I (M\^jQj\ v ) t 

+ 9 o _ r o c L<M|Ho) ; (3-79) 

where |H ) is defined by eq. ( |3.43|) and \v) stands for the valence quark state. 

The Hamiltonian in terms of the collective variables can be obtained from 
the Lagrangian eq. (|3.73|) as follows. The canonical momenta U a (a = 1, ■ • ■ , 8) 
conjugate to q a are 

dL 



n Q - — 

oq a 



- f d 3 r ^X a i) + ZmesonqA 

2 Jb 



+1' 'mesonQM^Ma + 2^^^ B ^ 8a ' (3.80) 

Here the quark field operator should be taken as the expectation value with 
respect to the rotated hedgehog ground state |H) of eq. ( |3.79|) by consistency 
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with the classical meson sector. Then, the canonical momenta are written as 

II a = (Xquark + T mC son) Qi^ia + quark + 2T meson) ?M^Ma H ^fi&ai (3.81) 

where the following expectation value for the quark field operator has been used 

^(H| J d 3 r ^\ a il> |H) = l quaxk qi5 ia + T ' quaik q M SMa + ^j-Y R 5 8a (3.82) 

where the following definitions apply 

T 3 _ |(m|A 3 |n)| 2 3 _ |(m|A 3 M 2 

X quark " 2^ P _ P + O ' 2^ r _ P ' 



3 x . |(m|A 4 |n)| 2 3 \(m\X 4 \n)\ 



2 



quark 2 ' ^ e° - uj° 2 ' ^ cj° - e° 
3 v |HA»| 2 

Z m,<>0 ^™ fc « 

= 1-B S . (3.83) 

We have summed over N c = 3 colors. By taking the expectation value of the 
quark operator with respect to |H), the classical Hamiltonian is obtained as 

H = (Mq Uark + M meson ) + -(Z quark + X mC son)9i 2 "^2^ quark meson )^ 

= Mo + | + § (3.84) 

with M° = M° uark + M^ cson , J = J quark + J mcson , and J' = X' quark + J' mcson . Quan- 
tization of the Hamiltonian can be done by promoting the canonical momenta to 
a quantum mechanical operator, 

# = M ° + 5G-?) J2 + 2^-^' < 3 - 85 > 

where C 2 , is the quadratic Casimir operator for flavor SU (3), J 2 the corresponding 
one for the spin of SU(2), and Yr the "right" hypercharge operator, needed to 
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represent the Wess-Zumino-Witten constraint, namely that physical states obey 



f^lphys) = |phys). (3.86) 

The eigenstates of the Hamiltonian can be written in terms of the Wigner— D 
functions |J0 



as 



= Jdim(p,q){a\D^(A)\b), 
\a) = \II 3 ;Y), 

\b) = \l'l' 3 ;Y>), (3.87) 

where (p,q) label the irreducible representation of SU(3), D^ p,q \A) the corre- 
sponding element, \a) and \b) the basis on which D(A) act, I the isospin of the 
baryon, J 3 the third component, Y the hypercharge and the primed quantities 
are the right isospin, right hypercharge etc. With this collective-coordinate wave 
function and |H) eq. ( p.79|) , the baryon is described by the wave function of the 
form 

|B) = ® |H). (3.88) 
The mass formula from eq. (|3.85|) and eq. Q5.51j) is 



M(p,q;II 3 Y:JJ 3 ) = M + i (± - 1 ) J(J + 1) 



+ ^(\{p 2 + pq + q 2 + S(p + q)}-j), (3.89) 



2PV3 ^ " " v " "J 4 
which yields the mass formulas for the baryon octet and decuplet 

M m = M» + ^ + ^r (3-90) 

Since the quark masses are ignored, all the particles in the baryon octet (decuplet) 
have the same mass. 
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3.6 The Gluons 



The treatment of the non-perturbative interaction between the pseudoscalar mesons 
and the quarks has been discussed in the previous sections. The chiral bag model 
contains besides the quarks and the octet pseudoscalar mesons other degrees of 
freedom, namely gluons and rf meson, which we will incorporate in a perturbative 
fashion. 

The gluons appear in two ways, as produced by the quark sources and through 
the vacuum properties of the cavity, i.e. the so called vacuum fluctuation. 

Let E a and B a be the color electric and magnetic fields, respectively. The 
index a denotes the gluon color and runs from 1 to 8. They satisfy generalized 
Maxwell equations for r < R, 

V-E a = J°' a , (3.91) 
V x E a = 0, (3.92) 

V B a = 0, (3.93) 
VxB" = J a , (3.94) 

with the boundary conditions due to confinement at the bag surface r = R, 

r-E a = 0, (3.95) 
r x B a = 0, (3.96) 

where J M ' a is the color charge current given by 

= ^7 M y^ (3-97) 

The boundary conditions resemble the case of the perfect conductor in electrody- 
namics, but the roles of E a and B a are interchanged due to the structure of the 
QCD vacuum ||. 
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The solution of Maxwell equations, for example, eqs. ( |3.91| ) and ( |3.92| ), can 
be written as 



= V / d 3 r G(r,r')J°' a (r'), (3.98) 
Jb 

with a proper static cavity propagator which satisfies the boundary condition |7(J 
[FT|j . Since all the valence quarks 4 have the same quantum numbers except for 
color, the color charge density operator, J^, and the current operator, J" al , can 
be written in the form 



4y = 9s<l>i(r)M r )^(<x\-f\PW a ai3, 

a/3 



J a val = ^t(r)a^(r)^(a|^|/5)at^, (3.99) 

a/3 " 

where \a) and \(3) (a,/5 = 1,---,N C ) denote the color states and (f) v (r) is the 
spatial, spin, and flavor wave function. 

With the help of eq. (|3.78j) , <fi v (r) becomes 

«r) = ^(r) + |- £ ™^ + f • £ ^ £ ^£. P.100) 

up to the lowest order in the collective variables, q a . Here <p° v h is the unper- 
turbed hedgehog quark state with K p = + and = in the lowest energy 
level. Because of the matrix element (n\Xi\v) ((m\\M\v)), the summation over 
n (m) is restricted to the hedgehog quark states with K p = 1 + and Mk = 0, ±1 
(strange quark states with K p = j p = | + and rrij = ±|). The conditions in the 
summation are necessary to be consistent with the particle-hole picture. 
Substituting the explicit wave functions to (f>® h , and 0^ leads to 



■#00 = ^EHyl/^W 

a,/3 



4 In case of — tt < 8 S < — 5, there is no valence quark. In this case, the quarks which occupy 
the lowest energy K p = + state are taken as valence quarks even though they are in the 
negative energy sea. 
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J? a i(r) = -9s^- (*x S)^E(«lf \®4flf>, (3-101) 

a,f3 

where S is the spin for K = defined by 

S, = -Xqt (3.102) 
in terms of collective variables ^ (i = 1,2,3) and of the moment of inertia X given 



in eq. (|3.84|) . The quantities p'(r) and /x'( r ) are given by 



p'( r ) = iVV 2 (j 2 (£^+ii(^)), (3-103) 



3I nk>o £ n- £ ° 1 Jo(e2r)ji(e°r) + ji(e°r)j (e°r) ) 

where the sum in fjf(r) runs over all positive energy eigenstates of K p = 1 + , e° 
is the energy of the n-th eigenstate, and 

fl n = N(3N 2 [ R r 2 dr( j (e v r)j (e n r) + 3l (e v r)j 2 (e° n r) ). (3.104) 



o 

Here N, aNi and (3N 2 are the normalization constants for the wave functions as 
appeared in eq. (|3.28| ); 

N- 2 = f R r 2 dr( j*(e° v r) + j*(e° v r) ), 
Jo 

N{ 2 = [ R r 2 dr(3l(e n r)+jl(ey)), 
Jo 

N 2 2 = [ R r 2 dr(f (e n r)+j!(e° n r)). (3.105) 
Jo 

Note that the color charge current can have a non-vanishing value due to the 
perturbation in the valence quark by the collective rotation. The color charge 
current does not get any contribution from the strange quark. 

Following the refs. @ and ff% , the color electric and magnetic fields have the 



form 



K,i = 0.^*E<«lvl0> fl i a * (3-106) 



a,/3 
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B 



val 



^(2M(r) + ^ R) Mr) 



47T 



with 



R 3 

p{r) 
fi{r) 

M(r) 



,t. 



47T 



R 



dr' (/(?), 



,,/3 



(3.107) 



These fields are of the same form as those of the MIT bag model || except for 
the numerical details due to the modification of the valence quark wave function 
by the chiral boundary condition. 



It is well known that the color electric field given by eq. (|3.107|) does not 
satisfy the boundary condition eq. (|3.95|) dynamically. One can have it satisfied 
by imposing that hadrons are color singlet states only at the level of expectation 
values ||. However, it is quite unnatural that while the color magnetic field in 
eq. ( |3.107| ) and all the other multi-pole electric fields automatically satisfy the 
boundary condition f72| , the monopole part requires an additional prescription. 
Here an alternative choice is proposed to make the monopole electric field satisfy 
the boundary condition. Suppose that a sphere of radius £ < E around the 
origin is excluded so that the S function term associated with the equation for 
the electric field, eq. ( |3.91| ), is not present. Then, the most general solution for 
this field is given by eq. ( |3.107| ) where now 



p(r) 



dr'p\r r ), 



(3.108) 



with an arbitrary A. The confinement boundary condition, eq. (|3.95| ), can be 
satisfied if A = R [21|. Because of the singularity of the field at the origin which 
introduces an additional 5 source in eq. ( |3.91|) , this function is not a solution to 
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the initial problem. However, if one is willing to accept that the electric field is 
discontinuous, i.e., zero at the origin, and assumes this function away from the 
origin, then all dynamical requirements will be satisfied. This solution satisfies 
the boundary conditions at the price of relaxing the continuity of the electric field 
inside the cavity. We classify the solutions of the electric field as; (solution I) if 
A = and (solution II) if A = R. 

We now enter the description for the Casimir effect. The vacuum fluctuation of 
the abelianized gluon fields is described by the time dependent Maxwell equations 
without any sources 

V ■ E a = 0, V • B a = 0, 

<9B a <9E a 
V x E a = — , V x B a = -— (3.109) 

dt dt ( J 

and satisfy the MIT confinement boundary conditions eq. ( |3.95| ) and eq. ( |3.96| ). 
The classical eigenmodes of the abelianized gluons can be classified by the total 
spin quantum number (J, M) given by the vector sum of the orbital angular 
momentum L and the spin S, 

J = L + S, (3.110) 

and the radial quantum number n. There are two kinds of classical eigenmodes 
according to their relations between the parity and the total spin; (i) M-mode 
with the parity P = — (— 1) J and (ii) E-mode with the parity P = — (— 1) J+1 . 
Here the extra minus sign is due to the negative intrinsic parity of the gluon. 

It is convenient to introduce the vector potentials, G a , and choose the Coulomb 
gauge condition; 

G a = and V-G a = 0. (3.111) 

Then, the electric field and the magnetic field are obtained in terms of the vector 
potential through the relations 

dG a 



E 1 ' 



dt 
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B a = V x G a . (3.112) 
Omitting the color index, from Maxwell equations eq. (|3.109 ), the solutions be- 



come 



fi) M — modes : 



^n,J,M 



(ii) E - modes : Gf J M 



AT M j J (u; n r)Yj i j iM (r 
Me 



(3.113) 



J 



2J+1 



jj+i(uj n r)Yjj +lM (r) 



+ \l „ T + \ jj-i(u n r)Y J>J _ 1)M (r) 



_ , ,_, M> _, , (3.114) 

where Yj ) j ) Af(r) are the vector spherical harmonics of the total spin J carrying 
the orbital angular momentum I. The energy eigenvalues are determined by the 
MIT boundary conditions eq. (|3.95| ) and eq. (|3.96| ) as 



(i) M - modes : X n fj(X n ) + jj(X r , 



0. 



(ii) E — modes 



3j{X n ) = 0, 



(3.115) 
(3.116) 



where X n = uo n R. The normalization constants Mm,e will be specified below. 

The field operator G(r, t) is expanded in terms of the classical eigenmodes in 
the form of 



G(r,f) = £ a M G M (r)e— »' + a\ u} G\ v} {T)e*** , . 
M v 



(3.117) 



where {v} denotes the quantum number set (n, J, M, X = E or M). 

The normalization constants Mm,e are determined in such way that the free 
gluon Hamiltonian operator 



becomes 



H = - ( d 3 r(E E + B B) 

2 J B 



M 



(3.118) 



(3.119) 



54 



when eq. ( |3.117| ) is substituted into eq. ( |3.118| ). It leads to a normalization 



condition for the classical eigenmodes given by 



/ B ^ G M- G M=db* Mw - (3 - 120) 



Then the normalization constants are determined explicitly 



U M = {X n i? 2 [j2(X n )-Jj-i(X n )j J+1 (X„) ] j" 1 / 2 , (3.121) 
U E = { X n R 2 j 2 ^(X n ) Y 1 ' 2 . (3.122) 



3.7 The jf Meson 



The rf field is incorporated in the Lagrangian eq. ( |3.1| ) by allowing the U field 
to be £7(3) valued. Since the rf cannot have any topological structure, it satisfies 
the usual Klein-Gordon equation of motion 



Moreover, the rf field decouples from the pseudoscalar octet meson fields. How- 
ever, there is some secondary coupling between them via the quark-?/ interaction 
on the bag surface. Notice that the introduction of rf field modifies the quark 
boundary condition to 



which shows how the rf field can affect the quark fields directly and the hedgehog 
solution indirectly. Assuming its effect to be small, however, the possible modifi- 
cation of the hedgehog solution by the rf field will be not considered. As in the 
gluon case the rf field will be treated perturbatively. 

The boundary condition for rf field arises from the continuity of the flavor 
singlet axial current on the bag surface, 



(d 2 - V 2 + m 2 K = 0. 



(3.123) 




(3.124) 



<H| r • ^775^ |H> = <H| f • (2f n V V ') |H). 



(3.125) 
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Before collective coordinate quantization, the hedgehog solution cannot have a 
flavor singlet axial current, so that the 7/ field is identically zero. When the 
hedgehog solution is rotated by the collective rotation, the matrix element of 
the flavor singlet axial current is linear in the spin operator. Thus, in order to 
satisfy the boundary condition eq. ( |3.125|) , the 7/ field should be linear in the 
spin operator. One possible static solution of the Klein-Gordon equation with 
this constraint is 

V ' = CS- V(^) = -C(S • r) (I+^e-V', (3.126) 
where the constant C can be determined by the boundary condition eq. Q3.125 ). 



56 



Chapter 4 
Anomalies 



4.1 Preliminary Remarks 

Symmetries and their corresponding conservation laws play an important role in 
describing the fundamental forces of nature. However, it might turn out that 
a certain conservation law or symmetry, which is valid in the classical level, is 
violated at the quantum level. This phenomenon is known as the anomaly. If 
the symmetry so violated is a local gauge symmetry, then such an anomaly must 
be cancelled at the quantum level. This is the so-called anomaly cancellation for 
gauge theories. We will see below that this is relevant in our development. If the 
symmetry is however global, then the anomaly can and does manifest itself in 
observables. A well known example is the Ua(1) anomaly. To see a role of Ua(1) 
anomaly in QCD, let's consider the QCD Lagrangian with the three light quarks 



Mann matrices for the color structure. If we take the current quark masses 
to be equal, the QCD Lagrangian is invariant under the global transformation 



(4.1) 



where D„ = - ig s § G» 



is the covariant derivative and A"'s are the Gell- 
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exp(iA a # a /2) in flavor space. Besides this symmetry, when the current quarks are 
massless, the QCD Lagrangian is also invariant under the global axial transfor- 
mation exp(i7 5 A a # a /2) in flavor space. In other words, the QCD Lagrangian has 
the chiral symmetry SUl(N f ) x SUr(N f ) with N F flavors. Through the spon- 
taneously symmetry breaking mechanism, the real symmetry of QCD becomes 
SUv(Np) and a pseudoscalar octet of Goldstone bosons appear. 

There are two additional global symmetries in the massless QCD Lagrangian. 
One is the global U(l) symmetry which corresponds to the conservation of the 
baryon number. The other is a global axial transformation [/^(l). The Ua{^) 
symmetry requires parity doublets in the hadron spectrum which are never seen. 
Therefore, it should be broken and there be the accompanying Goldstone boson. 
The only known candidate for this particle is rj . However rf is too heavy to be 
regarded as the Goldstone boson. This is known as the £7a(1) problem. The 
resolution of this problem may be in the fact that the Ua{^) is not a physical 
symmetry, i.e. the £7a(1) symmetry is broken explicitly due to a quantum effect, 
the so called Ua{^) anomaly. In addition to the Ua(1) problem, without the 



existence of the Ua{X) anomaly, the process tt — > 27 cannot be understood JF3 



We have been discussing, in the previous chapters, two phase scenarios in 
which the theory is described differently in each phase. We next show that in 
them the realization of the symmetries is more complex than in conventional field 
theory. The presence of two phases generates two sorts of anomalies, one global 
and the other local. The global symmetry involved is the Uy{X) corresponding 
to the baryon number and the local one is the local QCD color anomaly. Both 
should be conserved in a realistic theory. We described above how the baryon 
number is conserved in the two-phase picture. When considered on its own, the 
bag boundary induces the baryon charge to leak from the interior which can be 
interpreted as an effect of an induced axial current on the surface which leads 
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to an anomaly in the UyiX) current. This baryon charge leaked from the bag 
interior is picked up by the hedgehog pion that lives outside of the bag in such a 
way that the total baryon charge is preserved |H| . 

As explained below, the bag boundary induces the color charge to leak out 
also. In contrast to the baryon charge case, there is no topological field outside to 
absorb the color charge accumulated on the surface, this charge must be cancelled 
by a boundary condition. We will find that this requires the presence of a surface 
term that violates the local color symmetry. This means that that the classical 
action violates the gauge invariance which is rectified only at the quantum level. 
The way anomaly figures in this case is opposite to the global case mentioned 
above where a symmetry which is manifest at the classical level gets broken at 
the quantum level. 

In this chapter, the first section is devoted to a derivation of the {7a(1) anomaly 
in QED, to understand how it is generated by using the Schwinger's model in 
(1+1) dimension, and to extension of the C/a(1) anomaly to the non-abelian 
anomaly. In the next section, the relation between the C/a(1) anomaly and the 
flavor singlet axial charge, introduced previously, is described. The discussion on 
the color anomaly follows and completes this chapter. 



4.2 Axial anomaly in QED 

There are many methods to obtain the Ua(1) anomaly ]73|]. Here, the method of 
taking the divergence of the axial vector current of QED in the position space is 
considered since this method will be used again in the discussion that follows. It 
is well known that operator products at the same space-time point are singular. 
So the axial vector current consisting of two fermion field operators 

jg(x) = ^(x)Yl^(x) (4.2) 
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may be singular. The regularization of this current leads to the anomaly flT4" |. The 



point splitting method may be used to regularize the current operator. In this 
regularization the operators are separated by a small vector e M in the following 
way: 

jg(x, e) = $(x + e/2) 7 ^(x - e/2) exp (ie £^ dy v A v {y)) (4-3) 
and the regularized axial vector current is defined as 

£ ■ (x) =kmjg{x,e). (4.4) 

° e— >U 

Here, the exponential of the gauge field is introduced to keep gauge invariance. e 

should be sent to zero after all the calculations have been performed. The Dirac 
equations 

{i $ — m + e J£)ty = 0, 

Tp(i@ + m-e4) = (4.5) 

yield the divergence of the point-split axial vector current 

= 2imP(x } e) - iejg(x, e)e" {d v A„{x) - d^A u (x)), (4.6) 

where the definition of the point-split pseudoscalar density 

P(x, e) = $(x + e/2) l5 ij(x - e/2) exp (ie jT^ dy"A u (y)^j (4.7) 

has been introduced. Naively taking the limit e — > would yield the classical 
(partial) conservation law, 

dfjg = 2irmpy 5 ip (4.8) 

leading to the exact conservation in the massless limit. However, this procedure 
is incorrect since the operator jg(x,e) is singular. Let's consider the vacuum 
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expectation value of the second term with non-vanishing e: 



= e v tr(^ 5 ^S(x - e/2,x + e/2)) exp (ie jT^ dy"A p (yf), (4.9) 
where the fermion propagator in the external field A p has been introduced as 

S(x - e/2, x + e/2) = (0|T^(x - e/2)ij)(x + e/2)|0). (4.10) 
The fermion propagator can be expanded in powers of A p 

S(x-e/2,x + e/2) = S (-e) + te [ d A z S (x-e/2-z) 4(z)S (z-x - e/2) + • • • , 



(4.11) 

where Sq(x) is the free fermion propagator. The first free term vanishes because of 
the properties of the Dirac trace. Since the degree of divergence in the expansion 
decreases, only the linear term in A p survives. In momentum space the fermion 
propagator has the following representation 

S{x - e/2, x + e/2) ~ -ie f {dq)e*** f (dp)e ip€ -^— 4(q)^—^ , (4.12) 

J J p — m p— $ — m 

where the abbreviation, (dp) = t^t, has been used. Substituting this represen- 
tation into eq. (|4.9|) gives 



e"{0\jg(x,e)\0) = -ie (dq)e-^A x (q) (dp)e v e ip 



tr (^- m2) ((p- g)2 - m2 ) exp ri_ e/2 d vMv)y I 4 -") 

Because of 75, the trace generates a linear divergent term, 

e u (0\jg(x,e)\0) = Aee^ J (dq)e- iqx q p A x (q) (4.14) 



P 



x+t/2 



£ u / (dp)e tpt — " Q ^exp [ie dy?A p (y)), 

yr — m z ){ (p — q) z — m z ) \ Jx-e/2 
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where the trace of 7 matrices, tr757 M 7 Q 7 A 7 /3 = —4ie^ aXl3 , has been performed 
with the convention 5 0123 = — e 012 3 = 1. For the limit e — ^ 0, the integral over p 
becomes 

(dp)e,P V - ™ 2 )( (P- l) 2 ~ ™? ) ~* / (<ip) F e "" 



4- [»)' 

ide a J p 
1 e„ 



4 



3tt 2 e 2 ' 



(4.15) 



with the help of the following integral 



e 



tpc 



The symmetrization 



leads to 



lim— = V (4.17) 

e^o e 2 4 



-iee I/ (0|^(z,e)|0)| e _ + o(^ A4 (x) - <VU*0) 



e 2 



. -^ aA %A A (x)F^(x)lim^ 

= -^s a ^F aP F^, (4.18) 

Collecting the results and using the fact that the pseudoscalar P(x, e) is regular, 
the divergence of the regularized axial vector current has the form of 

dtfLgfr) = 2imP(x) - -^- 2 e a ^F a(3 F, u . (4.19) 

This equation, which expresses the non-conservation of the axial vector current 
even for case of the massless fermion, is known as Adler-Bell-Jackiw anomaly. [] 

1 When we apply the point-split regularization to the vector current, the term e !y tr(7 M S(x — 
e/2,x+e/2)) appears in the divergence of the vector current. Although tr('y li S(x — e/2, x+e/2)) 



is quadratically divergent, due to the symmetric structure as in eq.(4.17), there remains only a 



logarithmic divergence so that the divergence of the vector current vanishes. 
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It was proved that this anomaly is correct to all orders in perturbation theory for 

qed [g. 

As we have seen in deriving the axial anomaly, the Dirac vacuum (or sea) 
plays a crucial role. To see how the axial anomaly is generated from the Dirac 
sea J76|, let's consider 2-dimensional QED, the Schwinger's model, for simplicity. 
The Schwinger's model is composed of one massless fermion coupled to an abelian 
gauge field. The corresponding Lagrangian reads 

£ = $(i$+4)il>-~F^F^. (4.20) 

We assume that the fermion has unit charge. In two dimensions, we choose the 
Dirac matrices as 

7° = o"2, 7 1 = 7 5 = 7°7 1 = o"3, (4.21) 

where the oVs denote the usual Pauli matrices. From the fact that there is no 
mass term in the Lagrangian, chirality is a good quantum number. Therefore, 
the two component Dirac spinor has the form of 

*-(?.)■ 

where ipL,R are the eigenstates of 75, i.e. 75^,^ = ±ipL,R- Classically, there are 
two conserved currents: 

f = t^V, £ = ^75^, (4.23) 

as in four dimensions. To see the quantum effect on these currents, let's assume 
that the system has a finite length L and satisfies the boundary conditions 

A fl (t,x = 0) = A fl (t,x = L), 
iP(t,x = 0) = -ij;(t,x = L). (4.24) 
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In addition to these assumptions, we choose the Coulomb gauge so that A can 
be neglected and A\ is independent of x. Then, with the above 7 matrices, the 
Dirac equation becomes 



d 



d_ 

dx 



-A, 



tp = 0. 



(4.25) 



According to the boundary conditions the fermion wave function can be expanded 
into the Fourier series 



-iE k t 



L 



exp 



• 27r , , ■ > 
^— \ k + - )x 



(4.26) 



which yields the following energy eigenvalue for the L— and R— fermion eigen- 

states 



El 



2ty 



k + 



E, 



R 



2tt 



k + 



+ A U 



(4.27) 



with k — 0, ±1, ±2, • • -. Each type of the fermion has an infinite tower of energy 
levels. For A\ = the energy levels for L— and R— fermions are degenerate. If A\ 
is not zero, the levels split; the energy of the L— levels increases whereas that of 
the R— levels decreases. For A\ = 2tc/L, the original level structure is reproduced 
exactly as it should be because of gauge invariance. 

Now suppose that the system is in the vacuum state in which all negative 
energy levels are filled up and all positive levels empty with Ax = 0. Increasing 
the value of A\ from to ^ produces a L— particle and a R— hole. This situation 
is shown in Fig. ETI. Because the electric charges of the particle and the hole are 
2 If A\ changes by the hnite value, from A\ = 0, the Wilson loop, 



exp 



dxAi(x) 



has the same value as that of A\ 
gauge transformation. 



0. Therefore, A\ 



2 TV 

L 



is equivalent to A± = under a 
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Figure 4.1: Shift of the energy levels of the fermion due to a change of A\. 



opposite, the total electric charge vanishes under the change of A\. 

Q(t) = Jdx f(t, x)=Q L + Q R = 0, (4.28) 

where Ql,r are defined by 

Ql,r = J dx $l,rJoiPl,r = J dx iPl,r^l,r- (4.29) 

Consequently the vector current is conserved. On the contrary, the axial charges 
can have a non-vanishing value according to its definition: 

Qs = Jdx j$(t, x) = Q L - Q R = 2. (4.30) 

3 In Fujikawa's method [j77f , a nontrivial Jacobian appears in the path integral measure of 
the fermion field when a chiral transformation is performed. Using the eigenstates of the Dirac 
equation as basis, we can construct the Jacobian which contains the term 

which is equal to the difference of the zero modes of each chirality once the volume integration 
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We can rewrite this expression as follows: 



AQ 5 = —A Ax (4.31) 

7T 



and per unit time 

AQ 5 _ L AAi 

At ~ 7T At ' 

Using its definition, we obtain the relation for the axial current 



(4.32) 



d f 5 = -doA x (4.33) 

7T 

and finally we can write the anomaly for Schwinger's model in a Lorentz invariant 
form 

= -e^A v = ±-e^\ (4.34) 

For the explicit calculation that we will perform in the next chapter, we con- 
sider the flavor singlet axial current in QCD. It has the same anomaly equation 
except for the appropriate group theoretic factor. Applying the same regular- 
ization as in the preceding discussion and with gluon fields, = 4fGv, where 
4j-'s are the Gell-Mann matrices, the flavor singlet axial current, A ^, satisfies the 
following anomaly equation: Q 

gv A o = _^£ TrG QT 

E <";j-r. (4.35) 



Aix 



is performed. The Atiyah-Singer's index theorem gives the following result: 
J2 f d 4 x ^(x) 75 ^„(x) = ^ fd 4 x F, V F^ 

n,zero 

in the QED of the (3+1) dimensions. 

From now on, the notation, for the gauge field strength will be used for non-abelian 
case. 
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where the trace of the Gell-Mann matrices, Tr^|4rj = if> has been used in the 
last line and Np is the flavor number of quark. 

4.3 Application of the Ua(X) anomaly to the fla- 
vor singlet axial charge 

As discussed in Chapter 2, the flavor singlet axial anomaly supplies the key to 
understand the proton spin problem. In this section we give the relation between 
the flavor singlet axial charge and the singlet axial anomaly. From eq. ( |4.35| ) and 
the fact that G a G^ can be written in terms of the Chern-Simons current, 
as 

l -G%G^ = d,K\ (4.36) 

with the definition of 

K , = e ,upa G a ^ _ | /a6oG jG^ j (4.37) 

where f a bc is the structure constant of QCD, we can define the following conserved 
axial vector current for massless quarks (in the chiral limit): 

Al = 4-N F ^-K,. (4.38) 

In the gauge, f\ Gq(x) = 0, the gluon spin spin operator, S 9 , becomes 

S? = -e ljk Gid°G k a . (4.39) 

Moreover, in this gauge, the cubic term vanishes for the spatial components of 
K^. Thus, one finds the relation, 

Ki = -Sl, (4.40) 

5 The gluon spin and orbital angular are not separately gauge invariant and hence a choice 
of gauge is necessary. More detail is given in the next chapter 
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and its forward proton matrix element, 



(p,S\K»\p,S) = -S' t Ag, (4.41) 

where is the proton spin and Ag is the net helicity of the gluon along the 
proton spin. 

Due to the conservation of A ^ in the chiral limit, its forward proton matrix 
elements are independent of the renormalization scale and the form factor, a , 
defined through the matrix element 

<p,S\Al\p,S) = a°Sp (4.42) 

should correspond with the value in the quark-model, i.e., a = AS. From 
eq. Q4.38D and eq. Q4.41J) , therefore, one has the scale dependent flavor singlet 



axial charge in terms of AE and Ag as: 

a°(Q 2 ) = AJ:-N F -^Ag(Q 2 ), (4.43) 

as given in Chapter 2. 

Finally, let's consider the gauge dependence of the Chern-Simons current. For 
an abelian case like QED, the gauge transformation, 

A M (x)->A M (x) + «VV(x), (4.44) 

induces the change in the Chern-Simon current 

K^x) - K^x) - ^[d,A(x)]e^ pa F^(x). (4.45) 

While the Chern-Simons current changes, its forward proton matrix element (or 
expectation value) does not since the expectation value of F pa vanishes due to 
the derivatives in the definition of the field strength. Thus, the flavor singlet axial 
charge is gauge invariant for the abelian case. 
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On the other hand, for the non-abelian case as in QCD, the situation is more 
subtle. Under 

G^x) -> U(x)G fl (x)U-\x) + -(d^Uix^U-^x), (4.46) 

9s 

with the definition, G^ = G a -f, the change of the Chern-Simons current becomes 

K, - K, + -e^d^G^U^d^U) 
9s 

+ ^- 2 e, vpu Ti{U-\d"U)U-\d a U)U-\<rU)}. (4.47) 
39; 

The second term is a total divergence so that does not contribute to the forward 
proton matrix element. Although the third term can also be shown to be a 
divergence [73], but it cannot be discarded because of the non-trivial topological 



structure |79[] of QCD. As a result, its forward proton matrix element is not 
gauge invariant. To avoid these problems in our discussion, we will treat gluons 
as abelianized fields in the next discussions. 



4.4 Color anomaly in the chiral bag model [ |S3|] 

In the previous chapter, the boundary conditions for the confined gluons in the 
chiral bag model have been given as 

r • E a = 0, r x B a = (4.48) 

with r the outward unit vector normal to the bag surface and a the color index. 
These conditions mean that the color electric fields, E a , point along the surface, 
while the color magnetic fields, B a , are orthogonal to the surface. These are the 
usual MIT bag boundary conditions. We will show that the leakage of the color 
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charge resulting from the 7/ field at the bag surface, as discussed at the beginning 
of this chapter, makes these boundary conditions change. As a result, the color 
electric field has a component normal to the bag surface and the color magnetic 
field a component along the surface. 

The boundary condition for the quarks due to rf field at the bag surface is 
given by 

(77 • r + e i75,?(/3) )^ = 0, (4.49) 

where 77 = rj'/f^ = 7///0 an d P is a point on the surface. At the classical level 
this boundary condition makes the color current of the quark, confined inside the 
bag, to have no leakage outside. As a result, the color charge of the bag is a 
constant in time, 

(0|Q a |0) = - Jj(3 (0\j a ■ v\0) 

= -0 a j[d/3(O|^7-*y^|O) = O, (4.50) 

where we have made use of the quasi-abelian approximation for simplicity and 
E = dB. Once the calculation is completed, we will extract the non-abelian 
structure by inspection. 

The quantum correction to Q a due to the 77 can be obtained by introducing 
the gauge invariant point-splitting regularization to the color current operator 
as before. If we use a point-splitting in time direction and choose the temporal 
gauge condition, = 0, the regularized color current operator in the quasi- 
abelian approximation at the surface becomes 

j? eg (/3) = g${p + e/2) 7 ^ exp (ig s j** dfG^zfytf - e/2) 

= 0^(0 + C /2) 7 ^V(/J- e/2). (4.51) 
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Then, we have 



Q a = (0|Q a |0> = -g s lim / d/3 <0| ^(/3 + e/2) 7 • f-^(/3 - e/2) |0) (4.52) 



as the regularized expression. Using the boundary condition, eq. (f4.49|) , the effect 
of the i] appears explicitly as 



'? 

= lime / d/3 (0| + e/2) 7 ■ r^sW - e/2) |0>t)(/3), (4.53) 

where the relation 

e -ij5V(P+e/2) e i l5V (f3-e/2) = e -i 75 ( r? (/3+ e /2)-^(/3-e/2)) g -| fo03+e/2), 7?(/3-e/2)] 

~ 1 - ij 5 f]([3)e (4.54) 

up to 0(e), and point splitting in the time direction, have been used. Rewriting 
the vacuum expectation values in terms of the components gives 

Q a = -z| lime jf d/3 Tr ( 75 7 • ? y ^ + (/3 - e/2; /? + e/2)^(/3), (4.55) 

where S + ((3,f3') is defined by 

S+((3,(3')= lim S(x,x') (4.56) 

from the confined fermion propagator S(x,x'). The multiple reflection expansion 
method |80| produces the relation 



Q a = g s lim / d/3 (01 M + e/2) 7 • r—e~ l ^ +e/2) e l ^- e/2) ^((3 - e/2) 10) 

e^OJs 2 



S+(/3 - e/2; + e/2) = 1(1 + 1e ■ rU 5 )S(/3 - e/2; (3 + e/2) (3 - U, 1e ■ r) (4.57) 

in Euclidean space. Q Here t/g = e* 75 ' 7 . Assuming that 77 on the bag surface does 
6 The 7 matrices of Euclidean space are related to those of Minkowski space according to the 
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not depend on the location, but only on time, 7 we have 

= | lime £ d(3 E Tr( 75 y (l E ■ * + U 5 )S{(3 - e/2; (3 + e/2)^ 
= p s lime £ Tr( 75 y7 s -rS(/?- e/2; /3 + e/2)). (4.58) 



Here the eq. QQ7 ) and the boundary condition of the confined propagator in the 



footnote have been used in order to get the last result. Because of the structure 
of the 7 matrices in eq. ( |4.58|) the singular contribution arises from the first 



order term in the gluon interaction in the perturbative expansion of the confined 
fermion propagator. In the Minkowski space the confined fermion propagator is 

S(/3-e/2;/3 + e/2)~ / d 4 x S ((3 - e/2 - x){-ig sl ■ G b (x)\ b /2)S (x - (3 - e/2), 

Jb 

(4.59) 

where So(x,x f ) is the free fermion propagator. Substituting this result into the 
eq. (|4.58| ) leads to the integral in momentum space, 

dQa = -^lime£rf/3Tr( 75 7-ry5(/3-e/2;/3 + e/2) S 



di] 

= gl Km e jf dp Tr ( 757 - t~ f (dq) G b a (q)e^ + ^ 



rules: x# = x, x^ — ixq,Je = ~*7i74 = 7o> an d 75 = —1e1e1e1 A so that all 7 M become the 
hermitian matrices. The boundary conditions for the quark due to 77 becomes 

7 E • ii/j = and ipj E ■ r = -^U 5 . 

From these, the confined fermion propagator satisfies the boundary conditions 

U5(a)S(a, x ) = 7^ • rS(a, x ) and S(x, a)U§{ot) = — S{x, a)j E ■ f 

for any vector a on the boundary. 

7 One can generalize the result to a variation with a space-time dependent rj(x, t) by replacing 
the ordinary derivative by a functional derivative. 
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/ 



{dp)\l a -^—,e ip ' e )- (4.60) 



From the relations 



\a \b\ xab 

Tr( 7 5 7 i 7V7 CT ) = -4ie i(3acT , Tr( — — ) = — (4.61) 



2 2/ 2 

and the fact that for the limit e — > the integral over p becomes 



/ 



we arrive at the final result in the form of a surface integral 



dQ a g. 



2 



e°' QCT i d(3 Tid a G a M. (4.63) 



drj 4-7T 2 /s 

We can write this result in terms of the color magnetic fields in the quasi-abelian 
case asi 

d ^ = N F ^-J d(3B a -r, (4.64) 
drj Att z Job 

when there are iV^ massless quarks. Note that eq. ( |4.64p has been obtained with 
unconstraint radiated gluon. Since we consider a confine theory, there should 
be the additional factor of 1/2 in the final result by making gluon confined only 
inside the bag. Then, the result becomes 



dQ ° N F ' J 



2 



I dp& a -i. (4.65) 

JdB 



drj 8ir 2 Job 

Eq. (|4.65 ) means that color charge disappears from the bag due to the time 
variation of the 77 field which is the artifact of the effective approach as ours. 
In other words, the color charge confined at the classical level leaks out at the 
quantum level. A simple way to remedy this problem is to introduce a surface 
term of the following type as a counter term in the action to remove the artifact: 

2 

Sct~-N f ^£ dpG a B a -i v (P). (4.66) 

57T^ JdB 
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Chiral invariance, covariance, and general gauge transformation properties allow 
us to rewrite eq. ( 4.66|) in a general form 

2 

Sci = i^r- 2 i d/3 JCnpTriMJl -MJ), (4.67) 

0^7T JdB 

where is the Chern-Simons current defined by 



K" = r""" ■■(GlGlp - l -g s f abc G a u G b a Gl) ( 1.08! 



and is the outward unit four vector normal to the bag surface. Here, U = 
e *v /fo e m/f as given in the previous chapter. This counter term describes that 
the colorless rj field outside the bag interacts with gluons at the bag surface by 
the help of quarks. We expect that this term may be generated naturally in 
the effective action of the chiral bag model as the interaction term between the 
neutral pion and photons in the effective action of QED |81[ . 

Note that eq. ( f4.67[ ) is not gauge invariant on the bag surface because of 
the Chern-Simons current. Thus at the classical level, the Lagrangian is not 
gauge invariant. However, at the quantum level, the invariance is restored by 
the cancellation between the anomalous term, eq. ( f4.65|) , and the surface term, 
eq. flEBTp . 



Including this color anomaly phenomenon, the action for the chiral bag given 
in the previous chapter is generalized to 

S = Sb + S§ + Sqb 
S B = J^xfaw-^TrG^G^, 

Ss = f /b ( Tld ^ dflU + T^ m ?' [ Tr ( lnf/t - lnf/ )] 2 ) + ■ ■ • + S wzw, 

S 9B = \l dwin^U^+i^K^lntf -\nU)\. (4.69) 
2 JdB 16n 2 
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From the presence of the Chern-Simons term in the surface action, the bound- 
ary conditions for gluon fields at classical level are affected. In place of the MIT 
conditions eq. ( |4.48| ), we have instead 

r-E a = -^|f-BV, 



8* 2 f< 



o 



r x B a = x EV- (4-70) 

The fact that the color electric field can have a radial component contrary to the 
MIT conditions plays an important role on the proton spin problem. Besides, 
the rf mass can be estimated by using these boundary conditions and the axial 
anomaly [[82 . 
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Chapter 5 

Flavor singlet axial charge in the 
chiral bag model 



In this chapter, we present the calculation of the Flavor Singlet Axial Charge 
(FSAC) in the chiral bag model scenario and demonstrate how the Cheshire 



Cat Principle [[LJ, [15] operates for the observables that are not topological as 
baryon charge discussed in a previous chapter. In order to do so we need a 
specific formulation of the model through its equations of motion and boundary 
conditions. We should stress that although we are truncating the model, the 
model itself is quite general and represents low-energy dynamics of QCD. 

We recall that the equations of motion have been shown in Chapter 3 and the 
color boundary conditions have been introduced in Chapter 4. Our calculation 
will be carried out in the static spherical cavity approximation, that is, our bag, 
polarized along z— direction, is a static sphere of radius R dividing two regions 
of space in which the theory is implemented by QCD for r < R, and by an 
effective meson theory for r > R. Besides, we treat, throughout the calculation, 
the gluons as abelian fields, i.e., we work in the lowest order approximation to 
QCD in perturbation theory. 
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First we calculate the various static contributions to this observable, and then, 
after considering the gluon spin, we proceed to include the Casimir contribution 
which arises due to the change in the boundary conditions of gluons associated 
to the color anomaly. 

5.1 The formalism 

To obtain the FSAC, we need to calculate the matrix elements of the flavor 
singlet axial current. According to the Cheshire Cat Principle (CCP) discussed 
in Chapter 3, the flavor singlet current appears in the chiral bag model as the 
sum of two terms, one coming from the interior of the bag and the other from 
the outside, populated among others by the 7/— meson [] 

< = + (5-1) 

Here, 9b and are defined, as before, as 6# = 8(R — r) and = 1 — 
with R being the radius of the bag. We demand that the Ua(X) anomaly is given 
in the model by 

d,Al = £ G%G^ a Q B + 2Um 2 r>V e B . (5.2) 

a 

Our task is to construct a FSAC in the chiral bag model that is gauge-invariant 
and consistent with this anomaly equation. Our basic assumption is that in 
the nonperturbative sector outside of the bag, the only relevant {7a (1) degree of 
freedom is the massive r\ field. This assumption allows us to write 

K,b = = 2/^17 (5.3) 

with its divergence given by 

f = V.rn^. (5.4) 

1 From now on we will omit the prime and write simply r\ for this meson. 
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The immediate question to ask is: what is the gauge-invariant and regularized 
A B such that the anomaly equation, eq. (|5.2|), is satisfied? To address this 
question, we rewrite the current, eq. (|5.1|), absorbing the theta functions, as 

such that 

^(4Uq+4U) = 2 fnm 2 vV Q B , (5.6) 
d» A l,B G = — E^'B^, (5.7) 

~ a 

where the relation G° VLV G lLV,a = — 4E a ■ B a has been used. The sub-indices Q 
and G signal that these currents are written in terms of quark and gluon fields 
respectively. In writing eq. (|5.6j), we have ignored the up and down quark masses. 
Going to the static situation, appropriate to our discussion, eq. ( |5.7| ) can be 
written as 

V-A° = ^5>°-B fl e fl . (5.8) 

" a 

We should stress that since we are dealing with an interacting theory, there 
is no unique way to separate the different contributions from the gluon, quark 
and rj components. In particular, the separation we adopt, (|5.6f) and (|5.7| ), is 
non-unique although the sum is free of ambiguities. We find, however, that this 
separation leads to a natural partition of the contributions in the framework of 
the bag description for confinement we use. 

5.2 The quark contribution 

The quark current is given by 

a Bq = (5-9) 
78 



where tp should be understood to be the bagged quark field, which means that the 
quark field is a cavity mode before turning on the quark-gluon interaction. The 
quark current contribution to the FSAC is given by 



1 B C 



(P| / d 3 r^7375^ |P> 

JB 



(5.10) 



with the proton state, |P), given by eq. ([3.88 



The calculation of this type of matrix elements in the chiral bag model is 
nontrivial due to the baryon charge leakage between the interior and the exterior 
through the Dirac sea. But we know how to do this in an unambiguous way. A 



complete account of such calculations can be found in pi, [3T| 3^]. The leakage 



produces an R dependence, as shown in Fig. |5.1|, which would otherwise not be 
there in the matrix element, eq. (|5.10| ). It is worth stressing that, as seen in the 



Fig. for zero radius, that is, in the pure skyrmion scenario for the proton this 
matrix element vanishes. The contribution grows as a function of R towards the 
pure MIT result that technically is reached for infinite radius. The result of this 
calculation was presented in refs. [EH], EQ]. No new ingredient has been added. 



5.3 The meson current 

We can get the rj field contribution in terms of the quark contribution. From the 
boundary condition eq. ( |3.125| ) for the r\ field, 

(P| r- (^77 5 V>) |P> = <P| v-{2UVr,) |P), at r = R, (5.11) 

and the fact that the bagged current, ^775^, is divergenceless for massless quarks, 
we have the following identity: 

(P| / d 2 s r ■ (r 3 2/ w Vr ? ) |P> = (P| / d 3 r^j 3 ^) |P>. (5.12) 

J dB JB 
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1.0 




R(fm) 

Figure 5.1: Various contributions to the flavor singlet axial current of the proton 
as a function of bag radius : (a) quark contribution a° B ; (b) rj' contribution a° 
and (c) the sum. 
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Substituting 77 field given by 

, = CS.v(— ), (5.13) 
where S is the proton spin operator, into the left hand side leads to 

(P| / d 2 sr- (r 3 2/ ff V77) |P) = C/^5 3 [2(l + y„) + $ e -*>. (5.14) 

Here 

y v = m v R, (5.15) 

Therefore C is determined by 

C = 7^- [2(1 +y v ) +y$]e v «(P\ I rf 3 r(^ 7 375^) |P). (5.16) 



Consequently, the 77 contribution becomes 



(P| / #r2/.V 3 ,, |P) = „,, ! + 'f' ; <P| / <Pr$ l3lb i, |P). (5.17) 
J b lyx. -\- y^) + jb 



Thus we have the result that the r] contribution outside is entirely given in terms 
of the quark contribution inside and the i] mass. 

In Fig. |5.1| we show the radial dependence of several contributions. We show 



the results of |21], |20|, which arise from the charge leakage mechanism, and follow 
the quark distribution. We show also the contribution of the rj just calculated, 
by taking the quark current matrix element also from Since the 77 field 

has no topological structure, its contribution vanishes in the skyrmion limit. Our 
calculation illustrates how the dynamics of the exterior can be mapped onto that 
of the interior by the boundary conditions. We may summarize the analysis 
of these two contributions by stating that no trace of the CCP is apparent in 



Fig. |5.1| . Thus if the CCP is to emerge, the only possibility one foresees, is that 
the gluons do the miracle! 
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5.4 The gluon current A% G 



Understanding the FSAC and its implications in the present framework involves 
crucially the role of the gluon contribution, in particular its static properties and 
vacuum fluctuations, i.e., the Casimir effects. The calculation of the Casimir 
effects in the next section constitutes the principal aim of this work. 

Since we have assigned the anomaly to the gluon fields, eq. (|5.7|) , the gluonic 
axial current has the form 

K,Bo = ^WA^(^ a - ±g s f abc A? b A Xc ). (5.18) 

Note that this expression is not locally gauge invariant. In fact there is no gauge 
invariant dimension-3 vector operator with the gauge field alone. 

With this gluonic axial current, the boundary condition, eq. ( |5.11|) , changes 

to 

(P| r ■ (A° Bq + A° Bg ) |P> = (P| r ■ (2^Vr?) |P>, at r = R. (5.19) 

Using this boundary condition and the anomaly, eq. (|5.8| ), the total FSAC of the 
proton can be constructed. Since the quark is a divergenceless field, we rewrite 
eq. (|5.8| ) in the form 

V ■ (A° Q + A° Bg ) = N F ^ ]T E a • B a Q B . (5.20) 

^ a 

Integrating this equation with respect to the bag volume after multiplying r, the 
matrix element for proton becomes 

(P| / d 3 r(A° Bn +A° ) IP) = -(PUVjr— / d 3 r VE a -B a r |P) 

JB It JB a 

+ (P| j dB <m-{A Bq +A BG )v |P). (5.21) 
Using the boundary condition, eq. ( |5.19| ), we see that the second term yields the 
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contribution from r] outside the bag. Therefore, we have 

(P\J B d 3 r(A BQ +A° BG ) + J_d 3 rA° v |P> = -(P\N F ^ J^r £E a • B a r |P) 

-<P| [_d 3 r(2f 7r m 2 r)V )r |P), (5.22) 

and, since there is no explicit coupling between gluons and the r\ field at the tree 
level in the model Lagrangian, the first term on the right-hand side of the above 
equation may be considered as the gluonic contribution to the FSAC, namely 

a o = ( P |_^i f d 3 rr3 VE a -B a |P>. (5.23) 

IX J B „ 



This corresponds to the second term in eq. (|4.43 ) in Chapter 4. 

Let us proceed to calculate the gluonic contribution. We begin by dividing 
the gluon current into two terms according to their origin 

Ab g — ^G,stat + ^G,vac- (5-24) 

The first term arises from the quark and r\ sources, while the second is associated 
with the properties of the vacuum of the model. One might worry that this 
contribution cannot be split in these two terms without double counting. That 
there is no cause for worry can be seen in several ways. Technically, it is easy 
to show it in terms of mode creation and annihilation operators. One can also 
show this intuitively by making the analogy to the condensate expansion in QCD 



33|], where the perturbative terms and the vacuum condensates enter additively 
to the lowest order. 

Let us first describe the static term which effectively accounts for the mixing 
of the light quarks with the E a • B a of the anomaly. The boundary conditions for 
the gluon field corresponds in our approaximation to the original MIT one ||. 
The source for it is the quark current that appears in the equations of motion 
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after performing a perturbative expansion in the QCD coupling constant, i.e., the 
quark color current 

gsh^y (5.25) 

where the ip fields represent the lowest cavity modes. In this lowest mode ap- 
proximation, the color electric and magnetic fields are given by 

E " = S,~p(r) (5.26) 

B^4{^(3r..r-.) + (^ + 2M W ).} (5.27) 
where p is related to the quark density p' as0 



p(r,T)= ^ dsp'(s) (5.28) 
and p, M to the vector current density 

rr 

p(r) = / dsp'(s) : 
Jo 



M{r) = / ds 



R /*'(*) 



s 3 



As considered in Chapter 3, the lower limit T is taken to be zero in the MIT bag 
model — in which case the boundary condition is satisfied only globally, that is, 
after averaging — and r = R in the so called monopole solution ~ m 
which case, the boundary condition is satisfied locally. 

Now we introduce the rj field. We perform the same calculation however with 
the color anomaly boundary conditions given in Chapter 4, 

Nfq 2 N f q 2 
r ■ E a = -^fr ■ B a rj{R), i x B a = x E a V (R). (5.29) 

8tt Jo 87T j 



2 Note that the quark density that appears here is associated with the color charge, not with 
the quark number (or rather the baryon charge) that leaks due to the hedgehog pion. 
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In the approximation of keeping the lowest non-trivial term, the boundary con- 
ditions become 

? • = -^r • B^(i2) (5.30) 



r x B a stat = -^fv x E a g r,(R). (5.31) 



N F g 2 s 

Here and are the lowest order fields ^ given by ( |5.26| ) and ( |5.27| ) 
and f](R) is the ?7-meson field at the boundary. The r] field is given by 



„ (r) = _g^ s . f i±jv: e -^ (5 . 32) 



where the coupling constant is determined from the surface conditions, eq. ( |5.16|) , 



the results of refs. [21], |22| and this expression for the rj. 

Note that the magnetic field is not affected by the new boundary conditions, 
since points into the radial direction. The effect on the electric field is just a 
change in the charge, i.e., 

p stat (r) = p(r,T) + Pv (R) (5.33) 

where 

The contribution to the FSAC arising from these fields is determined from 
the expectation value of the anomaly 

N F a s 



a G,stat 



<P|-— / d 3 r r 3 EJ at • B« JP>. (5.35) 

7T JB 



The result of this contribution appears in Fig. |0| , where we show the MIT so- 
lution, the monopole one and the correction associated to both due to the color 
couplingF]. One sees that including the rj contribution in p s tat{ r ) produces a non- 
negligible modification of the FSAC but does not modify the result qualitatively. 

3 We have also investigated electric fields of the form + Br)f, but the results do not 
change much with respect to the ones shown since the B term tends to be small. 
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Figure 5.2: Dependence of aP G stat on the choice of V and the boundary conditions 
as a function of bag radius : (a) with an MIT-like electric field without rj coupling, 
(b) with a monopole-like electric field without r\ coupling, (c) with an MIT-like 
electric field with r\ coupling, and (d) with a monopole-like electric field with rj 
coupling. 

What is most striking is the drastic difference between the effect of the MIT-like 
electric field and that of the monopole-like electric field: The former is totally in- 
compatible with the Cheshire Cat property whereas the latter remains consistent 
independently of whether or not the rj contribution is included in p s t a t- 

Before going to the A^vac, let's consider the gluon spin to see another observ- 
able where the monopole solution seems to be favored by experiment. 
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5.5 The gluon spin in the chiral bag model [gjj] 

In this section we study the gluon polarization contribution to the proton spin T, 
which is identical to Ag in previous chapters, at the quark model renormalization 
scale in the chiral bag model. It is evaluated, as we will show later in detail, 
by taking the expectation value of the forward matrix element of a local gluon 
operator in the axial gauge A + = in the light cone frame. We show that 
the confining boundary condition for the color electric field plays an important 
role in the outcome. When a solution satisfying the boundary condition for the 
color electric field, the so called monopole solution, which is not the conventional 
one, but which we favor, is used, T has a positive value for all bag radii and 
its magnitude is comparable to the quark spin polarization. This results in a 
significant reduction in the relative fraction of the proton spin carried by the 
quark spin, which is consistent with the small flavor singlet axial current measured 
in the EMC experiments. 

As we have discussed in Chapter 2, the EMC experiment revealed the 
surprising fact that less than 30% of the proton spin may be carried by the quark 
spin. This is at variance with what one expects from non-relativistic or relativistic 
constituent quark models. This discrepancy — so called "the proton spin crisis" — 



can be understood as an effect associated with the axial anomaly [28]. If we follow 



the argument |33J] given in Chapter 2, the experimentally measured quantity is 
not merely the quark spin polarization AS but rather the flavor singlet axial 
charge, to which the gluons contribute through the axial anomaly. Therefore, to 
understand the proton spin crisis, the sign of the gluon contribution is crucial. 

Although not directly observable, an equally interesting quantity related to 
the proton spin is the fraction of spin in the proton that is carried by the gluons. 
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In Ref. [Q, the gluon spin T is introduced as 

\ = ±E + L Q + T + L G , (5.36) 

where £q,g is the orbital angular momentum of the corresponding constituent 
and T is defined as the integral of the polarized gluon distribution in analogy 
to E. The spin of course is gauge-invariant but the individual components in 
eq. (|5.36| ) may not be. T can be expressed as a matrix element of products of the 



gluon vector potentials and field strengths in the nucleon rest frame and in the 
A + = gauge. When evaluated with the gluon fields responsible for the N — A 
mass splitting, T turns out to be negative, T ~ — O.labag, m the MIT bag model 
and even more so in the non-relativistic quark model. 

By contrast, there are several other calculations that give results with opposite 
sign. For example, the QCD sum rule calculation MB yields a positive value 



2T ~ 2.1 ± 1.0 at 1 GeV 2 . In Ref. [|S6"|, it is suggested that the negative T of 



Ref. [55] could be due to neglecting "self-angular momentum." The authors of 
show that when self-interaction contributions are included, one obtains a positive 
value T ~ +0.12 in the Isgur-Karl quark model at the scale ~ 0.25 GeV 2 . 
Once the gluons contribute a significant fraction to the proton spin, due to the 
normalization eq. ( |5.36| ), the relative fraction of the proton spin lodged in the 



quark spin changes. Thus, the positive gluon spin seems to be consistent with 
the EMC experiment. 

We address this issue in the chiral bag model and pay special attention to the 
confining boundary condition for the gluon fields. 

Let us start by briefly reviewing how the gluon spin operator was derived in 



Ref. [55], |87|. From the Lagrangian 



£ = _Itt( G""G /( „ ) (5.37) 
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with = \G a ^ v , one gets the gluon angular momentum tensor 

Milv \ = 2Tr ^x u G^ a G a x - x x G^ - [x v g^ x - x x g^)C. (5.38) 
Integrating by parts, we have 

M „ux = 2Tr^-G^(x^ A - x A ^)A a + G^A^ + G^A A 

+d a (x»G, a A x - x x G, a A») + iG^G^(xV A - xV))(5.39) 

with A M = ^-A^. It seems reasonable to interpret the first term as the gluon 
orbital angular momentum contribution and the second as that of the gluon spin, 
while recalling that this is a gauge dependent statement. We will not consider 
the fourth term hereafter, since it contributes only to boosts. In Ref. p5| , 57| , the 
third term is also dropped as is done in the open space field theory. When finite 
space is involved, as in the bag model, dropping this term requires that the gluon 
fields satisfy boundary conditions on the surface of the region, as we next show. 
Let us express the gluon angular momentum operator in terms of the Poynting 
vector, i.e., 

J G = 2Tr [ d 3 r[r x (E x B)]. (5.40) 
Now doing the partial integration for B = V x A, we have 

J k G = 2Tr j^tfVjE^r x V) fc A/ + (E x A) fc ) - J^d 2 r{r -E)(r x A) fe | . 

(5.41) 

The surface term is essential to make the whole angular momentum operator 
gauge-invariant, but the surface term only vanishes, if the electric field satisfies 
the boundary condition on the surface, 

r ■ E = 0. (5.42) 

This is just the MIT boundary condition for gluon confinement. However, the 
static electric field traditionally used || does not satisfy this condition. Instead 
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the color singlet nature of the hadron states is imposed to assure confinement 
globally. 



We next show that the negative T of Ref. |55j results if this procedure to 
confine color is imposed. To proceed, we choose the A + = gauge and write the 
gluon spin in a local form as 

T = (P,t \2TiJ^d 3 x(^E x A) 3 + B ± ■ A_l)|P,T), (5.43) 

where _L denotes the direction perpendicular to the proton spin polarization and 
the superscript + indicates the light cone coordinates defined as x = ^(x°±x 3 ). 
We shall evaluate this expression by incorporating the exchange of the static gluon 
fields between i— th and j—th quarks (i 7^ j) which are responsible for the N — A 
mass splitting in the bag model. 

As discussed in Chapter 3, in the chiral bag model, the static gluon fields are 
generated by the color charge and current distributions of the i— th valence quark 
given by |H| 

^ a (r) = £p(r)| , (5.44) 
J?(r) = ^3(rxS)^|, (5.45) 

where p(r) and //(r) are, respectively, the quark number and current densities 
determined by the valence quark wave functions given in Chapter 3. They are 
very similar in form to those of the MIT bag model. There is, however, an 
essential difference, namely, that the spin in the chiral bag model is given by 
the collective rotation of the whole system while in the MIT bag it is given by 
a individual contribution of each constituent, i.e., there is no index % in the spin 
operator in eq. ( |5.45| ). 

The charge and current densities yield the color electric and magnetic fields 
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as 



where 



Bf = £ {S (w(r) + 4£> - ^) + K (r • S)^l f , (5.47, 



/i(r) = / dr'n'(r'), (5.48) 



M(r) = j R dr'^^-. (5.49) 



The quantity Q(r) can be determined from Maxwell's equations. The most gen- 
eral solution can be written as 

Q A (r) =4tt f r dr'r' 2 p(r'), (5.50) 

with an arbitrary A. The standard procedure is to choose A = so that the electric 
field is regular at the origin. This has been the adopted convention in the early 
days |§. We will refer to this solution as Qo(r). However, Qo(r) does not satisfy 



the local boundary condition eq. ( 5.42 ), since it is normalized as Qq(R) = 1. In 
Ref. ||, the fact that hadrons are color singlet states, had to be imposed in order 
to justify the use of this solution. 

Another solution, namely monopole solution presented in previous discussion, 



is obtained by setting A = R |21j and we will look for its consequences here. This 
choice satisfies the local boundary condition but requires the relaxation of the 



continuity of the electric fields inside the bag. It has been shown in [^TJ , and will 
be shown here again, that these two solutions, Qo(r) and Qr(t), lead to dramatic 
differences for certain observables. 

By using the static Green functions and the Coulomb gauge condition, one 
can obtain time-independent scalar and vector potentials from the charge and 
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current densities, eqs. ( |5.44| ) and ( |5.45| ), 

•rw = ("I) 

UfW = |»H(Sxr)| (5.52) 

where 

From these, the appropriate scalar and vector potentials satisfying the A + = 
gauge condition can be constructed: 

Af(r) = «&f(r), 

A?(r) = U?(r)-V ['dC*Kx,y,0, (5.54) 

Jo 

where the direction of the proton polarization is taken as that of the z— axis. 
Finally, we obtain 

#ja=l V JV 

+ f d 2 sz-r(V(x) [ Z d(Ef(x,yX) 

JdV \ JO 

-*7f(x) J o Z d(Ef(x,y,()))\PA) 
4 -a s \ R rdr Q x (r)(h(R) - 2h{r)), (5.55) 



3 

where a s = g^/An. The numerical factor in front of the final formula comes from 
the fact that SaW^baryon = —8/3 for i ^ j so that 

EE( P .Tl^H|P,T) = -2, (5.56) 

and the integration over angles yields 1/3. It is different from 8/9 of the MIT 
bag model |5J|, which comes from the expectation value 

EE(P>tl^f |l P 't) = - 4 /3. (5.57) 
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It is interesting to note that, if we naively substitute the static gluon fields 
and U" of eqs. ( |5.51| ) and ( |5.52| ) satisfying the Coulomb gauge condition into 

the second term of eq. fl5.41| ), we get 

4 r R 

V = --a s / rdr Q x (r)h(r), (5.58) 
3 Jo 



which is the same expression [plj p8| that was used in the previous section to 
evaluate the anomalous gluon contribution to the flavor singlet axial current a 
with the extra factor (—Npa s /27r), i.e., a = S — (A^?a s /27r)r / A . On the other 
hand, in Ref . |$6] , the gluon spin T instead of T' is used for the anomaly correction 
term because the calculation is performed in the A + = gauge. 

If the gluons can contribute to the proton spin, then the collective coordinate 
quantization scheme of the chiral bag model has to be modified to incorporate 
their contribution. That is because there is a natural sum rule namely that the 
total proton spin must come out to be |, whatever the various contributions 
are. In the chiral bag model, where the mesonic degrees of freedom also play an 
important role, the proton spin is described by the following contributions 

\ = |S + Lq + r + Lq + Lm, (5.59) 

where Lm, the orbital angular momentum of the mesons, has to be added to 
eq. Q5.36p . The proton spin is generated by quantizing the collective rotation 
associated with the zero modes of the classical soliton solution of the model 
Lagrangian. To the collective rotation, each constituent responds with the corre- 
sponding moment of inertia. The moments of inertia of the quarks and mesons, 
Xq and Xm, have been extensively studied in the literature |31|j . Substitution of 
the color electric and magnetic fields, given by eqs. ( |5.46| ) and ( |5.47| ) respectively, 
into eq. ( p.40[ ) defines a new moment of inertia of the static gluon fields with 
respect to the collective rotation as 

(J G ) = -X g lu, (5.60) 
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where the expectation value is taken keeping only the exchange terms, and u is 
the classical angular velocity of the collective rotation. 

We show in Figs. |5.7| (a) and (b) the gluon moment of inertia evaluated by 
using the color electric fields with Qr{t) and Qoir). In the case of Qr(v), X g 
is positive for all bag radii and comparable in size to Xq, the quark moment of 
inertia. On the other hand, Qo(r) results in a negative X G . This "negative" 
moment of inertia may appear to be bizarre but it may not be a problem from 
the conceptual point of view. The X G defined by eq. ( |5.60| ) can be interpreted 
as the one-gluon exchange correction to the corresponding quantity of the quark 
phase, which is still positive anyway. The point is that the spin fractionizes in 
the same way as the moment of inertia does. This means that we have 



1q 




vx G 


+x M y 












+ 1mY 




Xm 





L G + T = — , ^ , ^ (5.6i; 

Lm ~ 2{x Q + x G + x M y 

Each fraction as a function of the bag radius is presented in the small boxes inside 
each figure. Note in the case of adopting Qr(t) that at the large bag limit the 
proton spin is equally carried by quarks and gluons somewhat like the momentum 
of the proton. The negative X G obtained with Qo{r), thus, yields a scenario where 
the gluons are anti-aligned with the proton spin. 

The dashed and dash-dotted curves in Figs, (a) and (b) show the values 
for T\ and T' x . For comparison, we draw ~E by a solid curve. Note that, because 
of the difference in X G , even |S is different according to which Q\ is used. Again, 
both To and r' are anti-aligned with the proton spin. Note of course that the 
negative r' is apparently at variance with the general belief that the anomaly is 
to cure the proton spin problem. 
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To conclude, we show in Figs. [5]9] (a) and (b) the flavor-singlet axial current 
including the Ua(1) anomaly given by 

«o = E - ^T' y (5.62) 

For simplicity, we neglect other contributions to studied in previous discussion. 
They show that the positive V is consistent with the small measured in the 
EMC experiments. The radius dependence of each component may be viewed 
as gauge dependence both in color gauge symmetry and in the "Cheshire Cat" 



gauge symmetry discussed by Damgaard, Nielsen and Sollacher [89 



5.6 The Casimir effect on the FSAC due to the 
color anomaly 

Finally we proceed to study the term Ac, vac-, which arises from the so called 
Casimir effect associated with the anomaly. The vacuum in the cavity and the 
perturbative vacuum in free space are different due to the geometry of the cavity. 
This difference might lead to observable consequences and it has been considered 



for many observables and also for the quarks in our calculation |T[], [3T], [32] , but 
never for the gluons. We proceed next to describe the Casimir effect for the 
FSAC. 

The calculation of Casimir effects is in general complex and is plagued by 
divergences, which have to be properly taken care of. In order to clarify these 
issues we consider first the Casimir energy pD[ , which is well studied, and shows 
the structure of divergences. 

In the canonical quantization formalism of field theories in infinite (free) space- 
time, vacuum energies are divergent. By the Wick's normal ordering procedure, 
which is based on the fact that the physical measured energy is the difference 
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between the state energy and the vacuum energy, which acts as the energy origin, 
these divergences disappear. However, when one considers a theory in a finite 
space region with a boundary, the vacuum is changed with respect to the free 
one, the change depending on the geometry of the region. The Casimir energy is 
defined as the difference between the vacuum energy in the presence of a boundary 



and that of free space pi 



For example, let us consider the case of the free massless scalar field 4> in the 
region between two plates normal to the z— axis. The plates are separated by 
the distance L with the left plate at z = 0. If one imposes Dirichlet (Neumann) 
boundary conditions, namely = {d z <p = 0) at the plates, one obtains from the 
Hamiltonian density, by direct summation of modes, the vacuum energy density 
as 



3A 4 A 3 



7T 



2 



£ °' =2^2^4^- 1440^' (5 - 63) 
where superscripts D(N) and the negative (plus) sign in the second term denote 
the Dirichlet (Neumann) boundary condition. Here the cut-off A has been intro- 
duced to treat the divergent summation. The first term is the divergent vacuum 
energy for infinite space-time corresponding to the limit A — > oo. Therefore, the 
Casimir energy has the form of 

^ = tf*-* = *5Sri55£J- (5 - 64) 

Although there is an ambiguity due to the first divergent term, the plates feel a 
well defined attractive force, which can be obtained from this energy density. On 
the other hand, if one uses dimensional regularization instead of the cut-off, the 
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vacuum energy density does not have any divergent terms and is given by, 



e °' N —T^Aj)^—uMs< (5 - 65) 

where the Riemann's zeta function, defined by ((a) = J2^=i n ~ a , has been used. 
Thus, in this regularization scheme, the Casimir energy can be obtained without 



1 7T 2 



any ambiguity With this result, it was pointed out [93 that when electromag- 
netic fields confined between plates are considered, which can be regarded as two 
scalar fields, one with Dirichlet and one with Neumann boundary conditions, the 
vacuum energy density becomes 

so that the Casimir energy has the same form as that obtained by dimensional 
regularization. However, this result is a special case. In general, it is known that 
divergences appear in the Casimir energy and depend on the regularization used 

0- 

Similarly, the Casimir effect on any physical observable, O, may be defined 
by the difference between the vacuum expectation value of O with a boundary 
and that without boundary, i.e., 

Ocas = (0|O|0) boundary - (0|O|0) free , (5.67) 

and has a similar structure of divergences to those of the Casimir energy. 

The quantity that we wish to calculate is the gluonic vacuum contribution to 
the FSAC of the proton. It can be done by evaluating the expectation value 

N F a s 



(0 B | - / d\ r 3 (E a ■ B a )\0 B ) (5.6* 

7T JV 



4 By fiat power divergences are "killed" in dimensional regularization. This is the power of di- 
mensional regularization for renormalizable field theories but one should be aware of that power 
divergences play a physical role in effective field theories where the standard renormalizability 
requirement is not applicable. |92| 
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where \0b) denotes the vacuum in the bag. Note that since the electric and 
magnetic fields are orthogonal to each other, there is no contribution from the 
MIT boundary conditions and free space. The standard way to evaluate this 
expectation value would be to expand the field operators in terms of the classical 
eigenmodes that satisfy the equations of motion and the color anomaly bound- 
ary conditions, eq. ( |5.29|) . Although well-defined, this approach is technically 
involved. We have not yet obtained any quantitative results to report. In here, 
we shall proceed in the opposite direction. Instead of arriving at the CCP as in 
the standard approach, we shall assume the CCP and evaluate the Casimir con- 
tribution with the expression that follows from the assumption. The idea goes as 
follows. 

The CCP states that at low energy, hadronic phenomena do not discriminate 
between QCD degrees of freedom (quarks and gluons) on the one hand and meson 
degrees of freedom (pions, rj,...) on the other, provided that all necessary quantum 
effects (e.g., quantum anomalies) are properly taken into account. If we consider 
the limit where the 77 excitation is a long wavelength oscillation of zero frequency, 
the CCP asserts that it does not matter whether we choose to describe the 77, in 
the interior of the infinitesimal bag, in terms of quarks and gluons or in terms of 
mesonic degrees of freedom. This statement, together with the color boundary 



conditions, leads to an extremely and useful local formula [82 



E« • B"^ -^f^G 2 , (5.69) 

where only the term up to the first order in 77 is retained in the right-hand 
side. Here we adapt this formula to the chiral bag model. This means that 
the couplings are to be understood as the average bag couplings and the gluon 
fields are to be expressed in the cavity vacuum through a mode expansion. In 
fact, by comparing the expression for the 77 mass derived in [32 using eq. (|5.69|) 
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with that obtained by Novikov et al |J6 in QCD sum-rule method, we note that 
the matrix element of G 2 in Q5.69| ) should be evaluated in the absence of light 
quarks. This means, in the bag model, the cavity vacuum. That the surface 
boundary condition can be interpreted as a local operator is a rather strong CCP 
assumption which while justifiable for small bag radius, can only be validated a 
posteriori by the consistency of the result. This procedure is the substitute to 
the condensates in the conventional discussion. 
Substituting eq. ( |5.69| ) into eq. ( |5.68| ) we obtain 

<0 B |- — / d 3 rr 3 (E a -B a )\0 B ) 

7T JV 

yi ^-(p\Ss\p)(0B\ I d 3 r±G 2 r 3 r 3 \0 B ) 



N F a s ^ 


( N F g 2 




{ 8tt 2 




( N F g 2 




{ 8tt 2 


xE 

n • 


IV 



fo 

!J{R) (P\S 3 \P)(N?^V 



' fo 

B n - E; • E n )r 3 f 3 , (5.70) 

where we have used that r\ has a structure like (S-r)y(R). Since we are interested 
only in the first order perturbation, the field operator can be expanded by using 
MIT bag eigenmodes (the zeroth order solution). Thus, the summation runs over 
all the classical MIT bag eigenmodes. The factor (iV 2 — 1) comes from the sum 
over the abelianized gluons. 

The next steps are the numerical calculations to evaluate the mode sum ap- 
pearing in eq. ( |5.70| ): (i) introduction of the heat kernel regularization factor to 
classify the divergences appearing in the sum and (ii) subtraction of the ultravi- 
olet divergences. 

As given in Chapter 3, the classical eigenmodes of the (abelianized) gluons 
confined in the MIT bag can be classified by the total spin quantum numbers 
(J, M) given by the vector sum of the orbital angular momentum L and the spin 
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s, 

J = L + S, (5.71) 

and the radial quantum number n. In the Coulomb gauge there are two kinds of 
classical eigenmodes according to the relations between the parity and the total 
spin J: 

tt = -(-1) j 

G^W 1 ") = ^MjjMY JtJtM (r), (5.72) 
vr=-(-l) J+1 



(i) M-modes 



(ii) E-modes 



J 



2J+1 jj+i(unr)Yj, J+hM (r) (5.73) 



+ V 2J+ \ j J - 1 ( Unr ) Yj > J - 1 > M ( i ') 

where Yj^,m is the vector spherical harmonics of the total spin J composed of 
the angular momentum £ and jg(x) is the spherical Bessel functions. The energy 
eigenvalues are determined to satisfy the MIT boundary conditions as 

(i) M-modes : 

X n j' J (X n )+jj(X n ) = 0, (5.74) 

(ii) E-modes : 

Jj(X n ) = 0. (5.75) 
From the results in Chapter 3, the normalization constants Nm,e are specified as: 

-1/2 



Mi 



M 



X n R 2 (j 2 j(X n )-jj^(X n )jj +l (X n )) 
N E = \x n R 2 j) ; (.V„) 



The first step is to calculate the matrix elements 

Q M = 2 / d 3 r(B* {u} • B M - E* {l/} • E {u} )x 3 x 3 . 

100 



(5.76) 
(5.77) 

(5.78) 



From eq. ( |5.72| ), we obtain 



E M (r) 
B M (r) 



(+iuj n )J\f M j J (uj n r)Yj j j M (r), 



(5.79) 



J 



2J+ 1 



j J+1 (uj n r)Yjj + i M (r) 



+ V 2J+ 1 l Jj ~ 1 ^ nr ^ >Yj ' J ~ 1,M ^^ 



, (5.80) 



for the M-modes and the similar equations with E and B being interchanged for 
the E-modes. 

We encounter in the calculation the following angular integrals 



J dQY* je>M ■ Yj^ M rl. 



(5.81) 

By using that x| = (4/3)^^/51^0 + 1/3 and the Wigner-Eckart theorem, we 
obtain 

J cMY* AM ■ Y JAM r 2 3 = CJ/ (J(J + 1) - 3M 2 ) + i (5.82) 

where cj^ is a constant that depends only on J and I. We have to perform the 
summation over M, which runs from — J to J, which cancels the contribution of 
the first term, therefore we can take effectively 1/3 as the result of the integral. 
Finally, we obtain the matrix elements for the M-modes as 

1 Jo Xn x 3 dx 



Qn 



(5.83) 



3 X*\fi(X n )-jj- X (X n )jj +x (X n )\ 

In the case of the E-mode, we obtain exactly the same formula except the minus 
sign in front of it. (Note that the formulas for the electric field and the magnetic 
field are interchanged.) 

We have found that the matrix elements for the E-mode vanish up to our 



numerical accuracy as shown in Fig. 5.3. Here, the solid line is the spherical 
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(5.84) 



Figure 5.3: jj(x) and I{x) as a function of x. 
Bessel function jj(x) and the dashed line is the integral 

= J y 3 dy jj(y) - 2 j - 2J ^ 1 1 ij-i(y) 

We see that the zeroes of /(x) and coincide, thus showing that Q^^^n) — 0. 
The analytic proof is given the appendix. 

In order to regularize the mode sum, we introduce a heat kernel factor exp(— rX n ) 

S(r) = £(2J+l)Qi A ?e-^, (5.85) 

n, J 

where we have carried out the trivial sum over M and the vanishing E-mode 
contribution is excluded. 

Fig. p21 shows the numerical results of the sum up to X maa; =100, 150, 200, 



250 for the 40 values of r from 0.0025 to 0.1 with the step 0.0025. We can see 
that below r < 0.06 the convergence is poor. However, it is enough to see the 
presence of an 1/t 2 divergence. If we fit the data above r > 0.06, we obtain 

ai06i_0£8i6 + 0478 _ 00285T 

T Z T 

Apart from a possible logarithmic divergence, there are quadratic and linear di- 
vergences as we set r equal to zero. We shall remove these divergences following 
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Figure 5.4: Diverging properties of S(r) as a function of the heat kernel regular- 
ization parameter r. All the magnetic modes up to uj n R{= X n )=100 (solid circle), 
150 (solid square), 200 (solid diamond) and 250 (solid triangle) are included in 
the sum. 
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a procedure commonly used in the Casimir problem A caveat on this pro- 
cedure will be highlighted in the next chapter. Now if we neglect the logarithmic 
divergence, the best way to get rid of the quadratic and linear divergences is to 
evaluate 

S(r) + 2rS'(T) + ^r 2 S"(r) = J2C2J+l)QnA^-^X n + 0.5T 2 X 2 n )e- TX ". (5.87) 

~" n,J 



Fig. |57| show the results on this quantity for 80 values of r ranging from 
0.0025 to 1. We see that no serious divergences appear anymore. By fitting the 
convergent data with the above expressions for r, we obtain for the finite part 
of the sum 0.0478, from the cubic function fit, and 0.0456, from the quadratic 
one. These results are comparable to the finite term of the above naive fitting 
procedure (|5.86|) , which yielded 0.0478. 

Once we have the numerical value on the mode sum, the gluon vacuum con- 
tribution to the FSAC can be evaluated simply as 

4,c as = t = X l X T^eV X (0 -° 478) ' (5 - 88) 

where y(R) is related to a° Bg as 



'1 + m v R) 
[2(1 + m v R) + {m^Rf^m^Rf 



V(R) = ~ To7i I - ^ ;TL^ 5S? aL- (5-89) 



We have used N F = N c = 3, a s = 2.2, f = ^/N F /2f v , ~ 122MeV and m v = 958 
MeV. Our numerical results are given in Fig. |5.6| . The quarkish component of 
the FSAC is given by the sum of the quark and r\ contribution, a° Bq + a° ri and 
the gluonic component by a? G stat + aP Gvac . Both increase individually as the bag 
radius R is increased but the sum remains small, < a^ otal < 0.3 for the whole 
range of radii. 
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0.2 




Figure 5.5: S(t) — 2rS'(r) + \t 2 S"{t) as a function of r. The finite term of S{t) 
is extracted by fitting these quantities to a cubic and quadratic curves. 
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Figure 5.6: Various contributions to the flavor singlet axial current of the proton 
as a function of bag radius and comparison with the experiment: (a) quark plus 
i] contribution (a° BQ + (b) the contribution of the static gluons due to quark 
source (a,c,stat)i ( c ) the ghion vacuum contribution (a% vac ), and (d) their sum 
( a totai)- The shaded area corresponds to the range admitted by experiments. 
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Figure 5.7: The moment of inertia associated with the collective rotation as a 
function of the bag radius and the proton spin fraction carried by each con- 
stituents. In the calculation, we have used (a) the "confined" color electric field 
with Qr(t) and (b) the conventional one with Qo(r). 
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Figure 5.8: The gluon spin T as a function of the bag radius, (a) and (b) are 
obtained with the color electric fields explained in Fig. |5~7 . 
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Figure 5.9: The flavor singlet axial current function of the bag radius, (a) 

and (b) are obtained with the color electric fields explained in Fig. |577 . 
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Chapter 6 



Discussion 



This thesis deals with the description of hadronic phenomena from the perspective 
of QCD and low energy hadronic effective theories. The underlying goal has been 
the analysis of the realization of the Cheshire Cat principle in a realistic 3+1 
dimensional scenario. We have reviewed in this thesis its formulation and have 
studied an extremely subtle case, namely the FSAC, an observable intimately 
related to the study of the spin of the proton, that lends support to the notion of 
the Cheshire Cat Principle in QCD governing strongly interacting systems. Before 
embarking onto the discussion of the results, let us stress that the consequences of 
quantum effects, through the chiral anomaly, the color anomaly and the Casimir 
phenomena have played a major role in the successful completion of this work. 

In the previous chapter, we have considered various contributions to the FSAC 
of the proton. As presented in Fig. |5.2| , the FSAC, which arises from the quarks, 
the rj' meson and the (MIT or monopole) static gluons, fails to fulfill the Cheshire 
Cat Principle (CCP). For all cases the FSAC vanishes in the R — > limit. The 
missing ingredient which is required to restore the principle is the Casimir effect 
that accounts for the fact that the modes are not free but strongly interacting 
via the boundary conditions of the cavity. The calculated Casimir contributions 
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lead to a nonzero value of the FSAC even for zero bag radius. 

Our numerical results are shown in Fig. |5.(j| . Standard MIT bag parameters 
were used for the calculation. The quarkish component of the FSAC is given by 
the sum of the quark and r\ contributions, a% Q + a°, and the gluonic component by 
oP G stat + a° G vac . Both contributions increase as the confinement size R is increased 
but their sum remains small, < a^ otal < 0.3, for the whole range of radii, giving a 
value consistent with the experiment, a exp = a°(oo) = O.IOIqjq. It is remarkable 
that a(R = 0) ~ a(R ~ 1.5 fm), while each component can differ widely for the 
two extreme radii. 

We have shown that the principal agent for the observed small FSAC in the 
proton, in the framework of the chiral bag model, is the CCP. It is the CCP that 
assures the cancellation between the different contributions: the quarkish and 
the gluonic. Note that the separation used by us is arbitrary and has no physical 
meaning, only the sum is physical. 

For a small bag radius, both components are small, so the small net FSAC is 
inevitable. This is consistent with the observation that in the limit as R — > we 
recover the skyrmion description, which gives vanishing FSAC at leading order, 
subject to small modification by matter fields at higher order. For a large bag 
radius - a limit that corresponds to the MIT bag model - both the quarkish and 
the gluonic contributions are separately large but they cancel each other. Our 
assertion is that this cancellation is mandated by the CCP. We should recall again 
that the separation between the quarkish component and the gluonic component 



adopted in eqs. (|5.6| ) and ( |5.7|) is entirely arbitrary although the sum is unique. 

Whether the different components by itself are large or small has no physical 
meaning. Only their sum does. Different separations would lead to different 
scenarios leading to the same small value. These different separations are analogs 
to gauge choices in gauge theories as suggested by some authors (see, e.g., fTS]). It 
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is tempting to speculate that in some limit, the FSAC is exactly zero and the the 
small nonzero value corresponds to a departure from this limit. Understanding 
this limit would allow a unique separation of the components. 

One of the principal results of this thesis is that it is possible to have a nonzero 
value for the FSAC at R = and it is of the same size as that for large R. 
The reason for this nonzero value is intimately connected with the CCP, since 
it is the finite part of the gluon mode sum which normalizes the value of this 
contribution at the origin. Moreover, the color boundary condition provides us 
with a decreasing rf field contribution which changes softly as a function of R. 
While the effect of the surface color anomaly term is generally small for all radii, 
the finite nonzero value of the FSAC for R = is assured by the surface boundary 
term. Thus the violation of the CCP observed in the previous calculations at 
R = I^TJ p2| is neatly eliminated by the color anomaly boundary condition. 
More importantly, the monopole structure of the color electric field previously 
proposed is found to be required for the sign that comes with the important static 
gluonic contribution from the quark source. We believe that this cancellation is 
a manifestation in the bag scenario of the recently discovered one for QCD f97| . 
The MIT configuration would strongly violate the CCP. We are thus led to the 
conclusion that the CCP requires the monopole configuration for the color electric 
field. Whether or not this configuration leaves undisturbed other - successful - 
phenomenology was discussed in . 



In calculating the gluonic Casimir effect, we made the ah initio assumption 
that the CCP holds, an assumption which is expected to be valid for small bag 
radius. We then extended it, in accordance with the CCP, to all bag radii. We can 
justify this only a posteriori by showing that the CCP assumption is consistent 
with what one obtains. Note, however, that the gluonic Casimir effect is most 
significant for small R where it is needed for the CCP and plays little role for large 
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R. Thus our assumption is validated. It would of course be more satisfying if one 
could obtain the CCP as an output of the formalism, not put in as an input. To 
this end, we need to solve appropriate equations of motion for gluons to the color 
boundary conditions Q4.70Q . Since the color anomaly boundary conditions are 
generated by the quantum effect of quark, it is necessary to consider an effective 
action for gluons which contains the quantum effect resulting from integrating out 
quark fields. Although gluons have been treated as abelianized (or Maxwell) fields 
in the previous chapter, we should give special care in the abelianization because 
of the nontrivial topological structure of the QCD p8| . Therefore, we need to 
construct an effective action for the bagged QCD and see there are additional 
corrections besides the Maxwell terms in the effective action. 

We should mention a caveat left unspecified in Chapter 5 in regularizing this 
Casimir contribution. Since aP G vac vanishes when rj field is removed, the so-called 
"vacuum contribution" is duly subtracted in what we have computed. However, 
we have also explicitly subtracted quadratic and linear divergences appearing 
from the mode sum by resorting to a procedure used in most of Casimir-type 
calculation |99| , which, as far as we know, is physically reasonable but has not yet 
been rigorously justified from the first principles []. The same caveat applies to our 
calculation as it does to others. The finite term we have obtained might therefore 
be subject to additional finite corrections by procedures in the renormalization. 

We should also stress that our result is at best qualitative. A better treatment 
(such as a more realistic gauge coupling constant running with the bag size, a 
more accurate calculation of A^f c , etc..) might modify our results quantitatively. 
Even so, we believe it to be quite robust that the overall FSAC is small, ~ 0.3 
and that it is more or less independent of the confinement size. 

1 It is found that the specific structure of the divergent terms depends on regularization 
schemes whereas that of the finite terms does not |93fl . 
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Appendix A 

Angular momentum basis of the 
wave functions for the strange 
and the hedgehog quarks 



In this appendix, we present explicit expressions of the angular momentum basis 
for the strange quark (s-quark) wave functions, \j,m.j) K , which are eigenstates of 
J 2 and J z , and those for the hedgehog quark (/i-quark) wave functions, \K, m^i, 
which are eigenstates of K 2 and K z . Here k — ±1 and i runs over 1, 2, 3, 4. 

The basis \j,rrij) K can be constructed by combining the orbital angular mo- 
mentum basis \l,m,j) and the spin basis \\,m s ) s , which are eigenstates of the 
orbital angular momentum operators L 2 and L z and the spin operators S 2 and 
S z , respectively. According to the angular momentum sum rule, there are two 
types. 

i) J 1 — I + \ an( i k — 1: 
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1^— ^) |i,m i - l/2>[l/2, -1/2), 



2j 



(a.i; 



ii) j = I — | and k = —1: 



3 ~ m j 



+ 



2j + 2 

j + 777 j + 1 \ a 



|j,m i -l/2)|l/2,l/2) s 
|j,777,--l/2)|l/2,-l/2) s . 



(A.2) 



2j + 2 

The conventions of Edmonds [1100(1 for the Clebsch-Gordan coefficients have been 
used. One can see that these two type of solutions have opposite parity for a 
given j. 



Using the following identities [101 



cos6 Y™ = 
sine e^Y^ 
sine e-**YT 



I — m + 1 I + m + 1 
21 + 1 2/ + 3 



'ym _|_ 



I — m I + 777 



2/ - 1 2/ + 1 



Z + m + 1 Z + m + 2 



. 2Z + 1 21 + 3 
Z - 777 + 1 l-m + 2 
21 + 1 2/ + 3 



r /+l + 



•ym-1 _|_ 



Z — 777 Z — 777 — 1 



2Z-1 2Z+1 

Z + 777 Z + m — 1 



2Z - 1 2Z + 1 



(A.3) 



where Y^" 1 is the spherical harmonics, one can show that the above states satisfy, 



(A.4) 



The basis for the hedgehog quark, \K, rriK)i, can be obtained by combining 
the total angular momentum J and the isospin I. In this case there are four types 
of states given by, 

i)i = l = l + h K = j-A 



\K, m K )i 



/ K -m K + l 
\ 2K + 2 



\j,m K - l/2> 1 [l/2, +l/2) t 
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+ ( K \™l 2 + 1 ) * I J, m K + 1/2) t j 1/2 , - l/2) t , (A.5) 
ii)i = 2 (j = Z-± A = j + ±Y 



A" 1 

+ ( ~^"* ) 2 b,m* + l/2)-i|l/2, -l/2) f , (A.6) 



in) i = 3 j = /-i K = 



2 / ' 



\K,m K ) z = + l y\j,m K - 1/2)^11/2, +1/2), 

+ f g 1 VlJ. ™* + V2)-ill/2, -l/2) t , (A.7) 



2A + 2 

iv)i = 4^ = Z + |, K = j-\ 



\K,m K U = (^^) 5 |j,m i ,-l/2) 1 |l/2,+l/2) t 

+ ( K ~ mK Y\j,m K + l/2) 1 \l/2,-l/2) t , (A.8) 



V 2A^ 

where the subscript t has been used to label states in the isospace. 

Since K = I for % = 1,2 and K — I — 1 for z = 3,4, one can check that 
\K,mjc)i and |A^, m^) 2 have parity (— 1) K , whereas |A, 772^)3 and |A, m^-) 4 have 
parity — (— 1)^. 

Using the identity eq. ( |A.4 ), one can show that the states \K, m K )i satisfy the 



relation 

a ■ r \K,m K )i = -\K,m K ) i+2 . (A.9) 
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Besides, applying the operator r-r to \K, mxji where r are the Pauli matrices 
in isospace, we have the following relations; 



2j-2K 

t ■ r \K, 771a-) i = I 1 m ^)3 - 2 o 7? I i — 1^' m ^)4, (A.10) 



2j — 2K 

t v \K,m K ) 2 = 2K+l \K,m K ) A -2 y o ^ ; - | A", 771^)3, (A.ll) 



2j — 2K 

r-r\K, m K ) 3 = 2K+1 \K, m K ) 1 - 2 v - ^ ; - |X, m K ) 2 , (A.12) 





+ 1) 


2K + 


1 




+ 1) 


2K + 


1 


y/K(K- 


+ 1) 


2K + 


1 




+ 1) 



2j-2K 

T-i \K, m K ) A = - \K, m K ) 2 - 2 — - - |X, to k )i, (A.13) 

where we have used the fact that K in m^-) 13 is j — \ and K in 777.^)2,4 is 
j + |. These relations have been used to get the energy levels for the hedgehog 
quarks. 
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Appendix B 

Proof that the sum over E-modes 
is zero 



We have shown that the contribution of E-modes to the Casimir effect of the 
FSAC vanishes numerically in Chapter 5. In this appendix, we show an analytical 
proof of this result ||102|| . Let's consider I(x) given in Chapter 5; 



I(x) 



dy y 3 



fAv) 



j 



2J + 



jjj+i(y) 



J + i 

2J+ 1 



and the boundary condition for the E-modes 



jj(Xn) 



0. 



(b.i; 



where X n = u n R. To prove I(X n ) = 0, it is convenient to use formulae of the 
Bessel functions rather than the spherical Bessel functions. From the definition 



3i\ x ) 



we can rewrite I(x) in terms of the Bessel functions 



7T f x 

l Jo 



Jj+i/ 2 (y) 



j 



2J+ 1 
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^,7+3/2(2/) 



J+]_ T 2 f \ 

2J + 1 Jj-i/2{y) 



(B.2) 



0. 



and the boundary condition 

Jj+l/2\X n 

Using the following formula by Schafheitlin [|103 





(B-3) 



{fi + 2) jT dx x^ +2 J 2 {x) = (// + 1) (v 2 



dx x^ 3 2 {x) 



(B.4) 



+- 



z^ 1 (zJl(z) - ^ + l)Mz)) 2 + z* 1 (z 2 -v 2 + -^+ l) 2 ) J 2 (z) 



we can arrange I(x) in terms of Bessel functions, its derivatives and simpler 
integrals with a lower power in y. 

In the next step, we use the recursion relations of the Bessel functions in order 
to write I(x) in terms of the Bessel functions with an index J + \- In addition, 
to get more useful expressions, we decompose I(x) into the three terms 



h(x) 
h(x) 



IT f x 

2 Jo dy y2j2j +V 2 ^ 

t / dy y 2 J\ 
1 Jo 



'0 
7T J 



2 2 J 
vr J + 1 



J+3/2 



(y) 



2 2j+uo d yy J j-v^- 

Using Schafheitlin's formula (|B.4|) and eq. (B.3| ), we have immediately 



h \X r , 



7T 1 



xlj ,2 {x n ) + j{j + r 



' .7+1/2 



d v Jj+i/2(y) 



(B.5) 
(B.6) 
(B.7) 

(B.8) 



4(2 jo 

The second expression has a more complicated form in terms of Bessel functions, 



i.e. 



h{X n ) 



J 



4 2TTT { [ + 3J + 2] Jo dy Jj ^ y) 
2" 

T 
4. 



+ 2 X nJj+3/2 (X n ) - -Xlj'j +3/2 (X n )Jj +3 / 2 (X n ) 



1 

+-x, 

2 



X„ 



J Z -3J 
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J 



J+3/2 



{x n 



From the recursion relation of the Bessel functions [|104||, one gets 



Jj+s/2{y) = -J' J+ i/ 2 (y) + - r+ y 2 ' Jj+i/2(y), 



y 



which yields by the boundary condition eq. (|B~ 



Jj+3/2(X n ) - -Jj + i/ 2 {Xn)- 



Similarly 



Jj+3/2(y) 



J + 3/2 



Jj+3/2{y) + J J+ i/2{y), 



which means 



J.J+3/2(X n ) — y ~ Jj+l/2(X n ) 



So finally the following form for I 2 {X^) is obtained 
h(X, 

4 / J + 1 l 

J j+i/2(y) 



+ l^XnJ + X n 



+ [J z + 3J + 2] / dy 
Jo 



[B.9) 



rl 7 l Jj+i/2(y)Jj + x /2 (y) 



y 



(J + l/2) s 



y 



J j+i/2(y) 



For J 3 an intermediate result after applying Schafheitlin's reduction formula ( |B.3| ) 

reads, 



^3 (X r , 



1 3 /2 1 2 / 1 2 

+ 2^n^J-l/2pQ - 2^n^J-l/2(^n)^-l/2(^n) + gX n J j_y 2 (X n 



+ \x n {X 2 n - J 2 + J)J^ 1/2 (X n )} 
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Using this time 



Jj-l/2(X n ) - Jj +1 / 2 (Xn) 



J -1/2 T , 



-1/2 

it can be rewritten in the form 
n J+l 



J 



J+l/2 



(*») 



h(X n ) 



4 2J + 1 
+ [J 2 -J] 



J' 



J+l/2 



\Xn) 



-x* - -X n J + l -x n 



(B.10) 



dy 



Jj+i/2(y) + 



2 J + 1 T/ 



J' 



J+l/2 



(y)J J+1/2 (y) 



(J+V2)^ 

+ ^ 2 j+i/2 ly; 



Before actually adding up all three pieces, it is better to have a further con- 
sideration. The integrals involve three different types of terms, namely Jj +1 / 2) 
Jj+1/2 Jj+i/2 an d J j+i/2- However, these are not independent from each other due 
to the differential equation they fulfill. The differential equation for the Bessel 
functions reads 



d 2 J J+1 / 2 (y) 1 dJ J+1/2 (y) 



dy 2 + y 
and therefore we have 

f'X n 



dy 



+ 1 



(J+l/2) 5 



/ dy 
Jo 



J% V 2{y) 



f X n d 2 

= -J o dy J J+1/2 (y)—J J+1/2 (y) 

fX n ( i 

= / i - J j+^(y) J j+i/2(y) + 

jo i y 



(J + l/2) 5 



J 



J+l/2 



(v) ■ 



Using this identity, I 2 (X n ) , I 3 (X n ) simplify to, 



h{X n ) = 



7T J 



4 2J+ 1 



Jj+l/2\^n) 



i(X n 



x x 3 

2 " 



1 



X U J + X,; 



(B.11) 



[■x n r 2 J 11 

+ [J 2 + 3J + 2]/ dy Jj +1/2 (y) ~ —Jj + i/2(y)J'j +1/ 2(y) 

l y j j 
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h{X r , 



n J + l 



12 



>(x n 



42J+1 



— Xt, — —Xr,J ~\- —Xr, 

2 n 2 2 



(B.12) 



J j+i/2(y) + „ + 2 Jj+i/2{y)Jj + i /2 (y) 



y 



+ [J 2 -J] dy 
Jo 

Collecting the results, we have 

J(X„) = -|2J(J + l)jf o dy -J J+1/2 (y)J' J+1/2 (y) 



*1 
2 



n t/2 / y 
J J+l/2\ yy r, 



In the final step, we use the recursion relation again to rewrite the derivative 
of the Bessel function in the first term of the equation above in terms of the 
Bessel function as, 



Jj+i/2(y) = \ (Jj-i/2{y) - Jj+3/2(y)) 



and the integral formula ||104 



dy l -J p {ay)J q {ay) = ay ^lMM^lz ^VjM _ WW , 
y p 2 — q 2 p + q 



(B.13) 



Then, the first term becomes 



dy -Jj+i/2(y)J'j+i/2(y) 



l 



-X n J' 



12 



(Xn), 



/o y ■ •"• !./(•/ • l; 

and can be cancelled by the second term. Therefore, we have the final result 



I(X n ) = 0. 



(B.14) 
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